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ABSTRACT 


Seventy-five  measurements  of  the  differential  elastic 

electron-deuteron  scattering  cross  section  were  made  for 

2  -2 

values  of  g   ranging  from  0.05  to  0.35  fm    at  the  NPS  LINAG 

with  electron  scattering  angles  of  60   to  120  .   The  charge 

radius  of  the  deuteron  was  determined  to  be  r„   =  1.9641  + 

Ed 
0.0074  fm.   Values  of  the  neutron  charge  form  factors  were 

calculated  using  Lomon-Feshbach  Models  1,  5,  and  15  and 

relativistic  corrections  to  the  deuteron  wave  functions. 

Lomon-Feshbach  Model  15  gave  -~   G^  (0)  =  0.0  200  +  0.00  58  fm2 

dq     n 

in  closest  agreement  with  the  value  of  ;r  G.^  (0)  determined 

dg     n 
by  thermal  neutron  work.   A  Rosenbluth  plot  of  the  data  at 

2 

q   =0.2  shows  good  agreement  with  theory. 
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I.   INTRODUCTION 

A.   OUTLINE  OF  THE  EXPERIMENT  AND  INTRODUCTION  OF  FORM  FACTORS 

The  primary  purpose  of  this  work  was  to  measure  the  charge 

2  2 

form  factor  of  the  neutron,  G   (q"),  at  low  values  of  q*  (the 

n 

square  of  the  Lorentz  invariant  four-momentum  transfer  from 

the  scattered  electron  to  the  neutron)  in  order  to  resolve 

differences  among  several  earlier  experiments  [1,2,3,4]. 

Secondly,  by  comparing  the  G   obtained  with  the  use  of 

n 
several  different  deuteron  wave  functions  to  the  G   pre- 

n 
dieted  by  accurate  thermal  neutron  experiments,  one  can  select 

the  S  state  and  D  state  combinations  of  the  deuteron  wave 

functions  that  are  in  best  agreement  with  experiment. 

Thirdly,  measurements  of  the  deuteron  cross  section  at 

2 
low  values  of  q  permit  the  rms  charge  radius  of  the  deuteron 

to  be  determined  accurately. 

These  three  objectives  were  achieved  by  measuring  the 
elastic  electron  -  deuteron  differential  scattering  cross 
section.   The  elastic  electron  -  proton  cross  section  was 
measured  in  order  to  normalize  the  measured  deuteron  cross 
section.   The  electron  beam  was  produced  by  the  Naval 
Postgraduate  School  linear  accelerator  (LIN AC) . 

Since  the  primary  purpose  was  the  measurement  of  the  charge 
form  factor,  it  is  useful  to  consider  briefly  here  what  a  form 
factor  is.   In  Section  II  a  more  detailed  explanation  will  be 
given.   The  cross  section/-^— ,  for  a  relativistic  electron 
scattering  from  a  nucleus  that  has  zero  spin  and  magnetic 

10 


moment  and  consists  only  of  a  point  .charge  has  been  calcu- 
lated by  Mott  [5].   The  experimental  cross  section  is  re- 
lated to  the  so-called  Mott  cross  section  by 

where  G  is  the  form  factor.   If  the  nucleus  has  no  structure, 

2 
i.e.,  it  is  a  point,  G  =  1  at  all  values  of  q  .   If  it  has 

2 
structure,  G  =  1  for  q   equal  to  zero  and  G  is  less  than  one 

2  2 

for  q   greater  than  zero.   A  reason  why  G (q  )  is  called  the 

form  factor  [3]  is  seen  when  the  cross  section  for  elastic 

scattering  of  an  electron  from  a  nucleus  of  charge  Ze  with 

finite  extent  is  written  [6]  as 

(*£  )    =(d£)     rfr(r)e^?d3r  f 

Vd.a  /Exp    Vdxt  /Mott  U      ze       J   '    {X    Z) 

where  P(r)  is  the  nuclear  charge  distribution  and  c^q  is 
the  three-momentum  transferred  to  the  nucleus.   The  term  in 
brackets  is  the  form  factor  G  and  is  the  Fourier  transform 
of  the  charge  distribution.   If  one  calculates  the  theo- 
retical cross  section  for  elastic  scattering,  it  is  seen 
that  the  charge  distribution  involves  the  product  of  the 
initial  wave  function  of  the  scatterer  at  rest  in  the  labo- 
ratory frame  and  the  final  wave  function  of  the  scatterer 
after  it  has  absorbed  the  recoil  momentum  c^q.   In  Section 
II  it  will  be  shown  that  a  relativistic  correction  to  the 
deuteron  wave  function  is  required  because  of  its  recoil 
velocity.   This  is  a  small  (a  few  percent)  but  important 
correction  in  this  work. 
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If  we  consider  elastic  scattering  of  an  electron  from 
an  infinitely  massive  (that  is,  a  nucleus  which  does  not 
recoil)  nucleus  of  charge  Ze  Equation  I--2  will  correctly 
describe  the  scattering;  P(r)  is  then  the  static  charge 
distribution.   However,  for  light  nuclei  where  recoil  is  im- 
portant, such  as  deuterium  and  hydrogen  which  were  used  as 
targets  in  this  work,  the  intuitive  concept  of  the  form 
factor  as  the  Fourier  transform  of  the  charge  distribution 
is  not  unambiguous.   For  example,  if  the  charge  distribution 
of  the  proton  were  spherical  in  one  frame  of  reference  then 
in  another  frame  the  Lorentz  contraction  of  the  distribution 
would  cause  it  to  be  nonspherical.   To  describe  form  factors 
in  a  way  more  useful  in  theoretical  work,  they  are  considered 

as  invariant  functions  of  the  four-momentum  q  rather  than  the 

"•  222     22  ~*2  2 

three -momentum  c^q,  where  ^  c  q   =  ft  c  q   -  (E.-Ef)   and  E. 

and  E.p  are  the  initial  and  final  energies.   When  the  recoil 

2  _  -2 
is  small  then  q   ^  q   and  the  concept  of  the  form  factor  as 

the  Fourier  transform  of  the  charge  distribution  has  useful- 
ness.  As  will  be  seen  in  Section  II  the  invariant  form 
factors  have  other  theoretical  meanings  in  terms  of  strong 
interaction  dynamics. 

There  is  another  complication.   Particles  and  nuclei  may 
have  magnetic  moment  distributions  as  well  as  charge  distri- 
butions.  These  distributions  may  or  may  not  have  the  same 
form  as  the  charge  distributions.   Thus,  a  particle  or  nuc- 
leus that  has  a  magnetic  moment  (but  no  quadrupole  or  higher 
moments)  will  have  to  be  described  by  two  independent  form 
factors.   The  decomposition  into  electric  and  magnetic 

12 


distributions  depends  upon  the  frame  of  reference.   For  ex- 
ample, the  elastic  scattering  of  an  electron  from  a  massive 
nucleus  with  charge  Ze  and  no  magnetic  moment  will  be  cor- 
rectly described  in  the  laboratory  frame  by  Equation  1-2  if 
the  nucleus  is  at  rest  in  the  laboratory  frame.   If,  however, 
the  calculation  is  performed  in  the  electron's  frame  of  refer- 
ence then  the  moving  nucleus  will  also  have  a  magnetic  field 
and  thus  an  apparent  magnetic  moment.   In  this  work  we  are 
ultimately  interested  in  the  electric  form  factor  of  the 
neutron.   A  particularly  useful  frame,  the  Breit  or  "brickwall 
frame",  separates  the  electromagnetic  interaction  of  the  nuc- 
leolus into  two  form  factors,  G^,  the  electric  form  factor, 
and  GM,  the  magnetic  form  factor.   Sachs  et  al  [7,8]  showed 
that  these  are  particularly  useful  form  factors  since  G   and 
G   can  be  expressed  as  the  Fourier  transforms  of  the  spatial 

distributions  of  charge  and  magnetization  in  the  Breit  frame. 

2    ~*2 

The  Breit  frame  is  the  frame  in  which  q  =  q   [7],  i.e.,  that 

has  E.  =  E.p  and  that  has  the  absolute  value  of  the  momentum 
of  the  incoming  nucleoli  the  same  as  the  absolute  value  of  the 
momentum  of  the  outgoing  nucleon  and  oppositely  directed  [9]. 

Consider  the  electric  form  factor  of  a  nucleon  in  the 
Breit  frame,  _  _ 

GE(q2)  =  $  p(r)  ±\±*   , 

o 

which  for  a  spherically  symmetric  particle  may  be  written  as 
GE(q2)  =  |~  $   sin(qr)r  p(r)dr. 
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2 

For  small  values  of  q   the  sine  may  be  expanded  to  give: 

GE(q2)  =_GE(0)  -  iq2  r2ras  ,  (1-3) 

where  r     is  the  root  mean  square  radius  of  the  charge 
distribution.   A  similar  expansion  may  be  made  for  the  mag- 
netic distribution  to  obtain  the  magnetic  r    .   The  proton 

rms       v 

and  neutron  each  have  different  form  factors.   Thus,  four 

form  factors  are  required  to  describe  them:   G   ,  G   ,  G   , 

P    n    p 
and  GM  .   The  r     obtained  by  using  Equation  1-3  is  deter- 

n 
mined  in  the  Breit  frame,  but  since  Equation  1-3  is  only  good 

2 

at  low  values  of  q  ,  the  r     value  determined  is  generally 

^         rms  r        J 

quoted  as  the  nucleoli's  charge  radius. 

B.   METHOD  OF  OBTAINING  THE  NEUTRON  ELECTRIC  FORM  FACTOR,  G 


EN 


The  process  of  obtaining  G    from  the  measured  elastic 

n 


electron-deuteron  cross  section  will  be  outlined  here  and 
discussed  in  detail  in  Section  II.   First,  the  elastic  elec- 
tron-proton cross  section  was  measured  at  a  given  scattering 
angle  and  scattered  energy  using  hydrogen  gas  as  a  target. 
The  hydrogen  was  then  replaced  with  deuterium  gas  and  the 
elastic  electron-deuteron  cross  section  was  measured  at  the 
same  scattering  angle  and  with  the  same  scattered  energy. 
The  incident  electron  energy  was  adjusted  to  give  this  con- 
dition.  All  geometries,  beam  currents,  gas  pressures,  and 
operating  procedures  were  kept  as  constant  as  possible  while 
both  the  proton  and  the  deuteron  data  were  taken.   In  fact, 
the  LINAC  was  operated  continuously  from  the  time  an  experi- 
ment started  until  it  was  completed  with  both  cross  sections 
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measured.   The  time  required  for  such  a  "run"  varied  from 
about  15  to  40  hours. 

Because  absolute  cross  sections  cannot  be  measured  here 
owing  to  systematic  uncertainties,  the  ratio  of  the  cross 
sections,  cr  /or ,    was  formed.   In  this  ratio,  many  systematic 
uncertainties  cancel  since  the.  cross  sections  are  measured 
under  the  same  conditions.   The  absolute  proton  cross  section 
has  been  measured  in  previous  experiments  [10,11,12],  and 
these  data  have  been  fit  by  deVries  [10]  to  formulas  from 
which  the  cross  section  can  be  calculated.   When  the  experi- 
mental ratio  above  is  multiplied  by  the  absolute  proton  cross 
section,  the  absolute  deuteron  cross  section  is  obtained. 

Several  steps  are  required  to  extract  G    from  this  cross 

n 
section.   The  experimental  cross  section  is  first  divided  by 

the  Mott  cross  section  of  the  deuteron  to  obtain  the  deuteron 
form  factor,  G. .   G,  accounts  for  all  the  structure  in  the 
deuteron,  i.e.,  the  electric  and  magnetic  distributions 
within  the  neutron  and  proton  as  well  as  the  convection  of 
charge  due  to  their  motion  within  the  deuteron  and  the  mag- 
netic moment  caused  by  this  convection  of  charge.   Following 
Jankus  [13] ,  the  deuteron  was  treated  as  though  its  two  con- 
stituent particles,  the  neutron  and  proton,  were  point  par- 
ticles.  The  distribution  of  charge  predicted  by  the  deuteron 

wave  function  is  then  contained  in  a  term  called  the  deuteron 

2 
charge  structure  factor,  which  depends  only  on  q   and  the  wave 

function  of  the  deuteron  ground  state.   E.L.  Lomon  [14]  has 

provided  three  selected  sets  of  values,  from  his  larger  group 

of  sets  of  values,  for  these  deuteron  charge  structure 

15 


factors.   Each  of  these  is  based  on  a  different  deuteron  wave; 
function  obtained  from  nucleon-nucleon  data  [15]. 

Drickey  and  Hand  [3]  experimentally  verified  a  proportion- 
ality between  the  magnetic  form  factor  of  the  deuteron  and  its 
charge  form  factor,  called  a  "scaling  law".  That  relationship 
was  used  to  remove  that  part  of  the  deuteron  form  factor 
caused  by  the  intrinsic  magnetic  moments  of  the  proton  and  the 
neutron  as  well  as  that  part  of  the  deuteron  magnetic  moment 
caused  by  the  convection  of  charge. 

At  this  point  only  the  extended  charge  distribution  of 
the  proton  has  not  been  accounted  for  in  the  extraction  pro- 
cess.  This  distribution  is  accounted  for  by  the  proton  elec- 
tric form  factor  derived  from  deVries'  data  (see  Appendix  A). 

The  neutron  electric  form  factor,  G„  ,  has  thus  been  extracted 

n 
from  the  deuteron  cross  section. 

The  above  procedure  is  justified  since  the  deuteron  is  a 
weakly  bound  structure,  i.e.,  the  proton  and  the  neutron  spend 
most  of  the  time  outside  the  range  of  nuclear  forces  and  one 
can  adequately  describe  the  charge  and  current  distribution 
within  each  nucleon  as  if  they  were  free  [13] . 

Early  in  this  work,  before  the  gas  target  system  was 

operating,  efforts  were  made  to  measure  G_,   by  a  process  simi- 

En 
lar  to  that  described  above  but  using  polyethylene  targets  of 

CH„  and  CD„  to  obtain  the  cross  sections.   There  were  several 
problems  with  that  work  but  the  most  important  limitation  was 
the  scattering  from  carbon  which  complicated  the  data  reduc- 
tion.  The  solid  target  data  which  were  obtained  were  not  as 

reproducible  as  the  gas  target  data  and  are  not  included  in 

this  work. 
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C.   RESULTS  OF  PREVIOUS  EXPERIMENTS  AND  OF  THIS  WORK 

In  this  work  the  neutron  charge  form  factor,  G   ,  was 

E 

2  n   -2 

measured  for  values  of  q   ranging  from  0.05  to  0.35  fm 

2 

The  slope  of  Gn   versus  q  was  obtained  by  a  linear  least 

En 
squares  fit  to  the  data. 

2 

The  slope  in  the  limit  that  q  _>  0  has  been  measured  in 

many  experiments  [16,17,10]  which  measured  the  interaction  of 
thermal  neutrons  with  atomic  electrons.   The  most  accurate 
and  most  recent  measurement  of  this  kind  is  that  of  Krohn 
and  Ringo  [4] ,  who  obtained 

-— -  Gp  (0)  =  0.0193  +  0.0004  fm2.  (1-4) 

dq    n 

This  is  probably  that  best  determined  of  all  nucleoli  form 

factor  parameters  [19]  .   Drickey  and  Hand  [3]  measured  G,., 

En 
by  measuring  the  ratio  of  the  elastic  electron-deuteron 

cross  section  to  the  elastic  electron-proton  cross  section 

using  liquid  hydrogen  and  deuterium  for  targets  and  found 

2  -2 

that  G^   =0.00+0.01  over  a  range  of  q   from  0.3  to  2.2  fm  . 
E  — 

n 
That  is,  they  found 

dq     n 
which  is  in  disagreement  with  the  thermal  neutron  measurements 

Casper  and  Gross  [20]  developed  relativistic  corrections 
to  electron-deuteron  scattering  and  reanalyzed  the  data  of 
Drickey  and  Hand.   Using  deuteron  wave  functions  derived  from 
the  Feshbach-Lomon  nucleon-nucleon  model  they  found  that  the 
data  of  Drickey  and  Hand  indicated 


~^T  GEr,(°)  =  0.0153  +  0.0067  fm2, 
dq     n 
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which  was  consistent  with  the  thermal  neutron  slope. 

Bumiller,  et  al.  [21,2]  extended  the  elastic  electron-deuteron 

2  2-2 

scattering  data  to  q   of  0.10  ^  q   ^  0  .  8  f  m   ,  using  Clip,  CD~, 

and  C  targets  and  substantiated  the  conclusions  of  Casper  and 

Gross  by  obtaining  a  slope  of  0,017  9  +  0.003G.   This  quoted 

error  of  +  0.0036  is  incorrect.   It  should  have  been  quoted  as 

+  0.0089.   Using  the  same  equipment  as  Bumiller,  Mader  [1] 

continued  the  elastic  electron-deuteron  scattering  work  but 

used  a  gas  target  to  avoid  the  carbon  subtraction  procedures 

which  were  necessary  with  the  solid  targets  used  by  Bumiller 

et  al.   Although  Mader' s  data  are  as  accurate  as  Refs.  21  and 

2,  his  measurements  showed  that  y  Gp  w  0  over  a  range  of 

dq     n 
2-2 
0.05  ^  q   ^  0.60  fm   .   In  this  work  we  have  resolved  this 

discrepancy  by  measuring  the  G„   slope  to  be 

n 

-~   G^  (0)  =  0.0200  +  0.0058  fm2, 
dq     n 

which  then  compares  favorably  to  the  value  of  0.0193  +  0.0004 

2 

fm   as  measured  directly  by  the  neutron-electron  scattering 

experiments. 
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II.   THEORY 

A.   ELASTIC  ELECTRON-PROTON  SCATTERING  FROM  A  POINT  PROTON 

Although  it  has  been  proved  [22]  that  the  proton  is  not  a 
point  particle,  it  is  convenient  initially  to  solve  for  the 
cross  section  for  elastic  electron-proton  scattering  from  a 
point  proton. 

The  initial  assumptions  are:  1)  the  proton  has  a  point 
charge  and  a  point  Dirac  magnetic  moment,  2)  the  one  virtual 
photon  exchange  adequately  describes  the  scattering  process, 
and  3)  the  electron  is  a  point  particle.   The  Feynraan  diagram 
describing  this  process  is 


where  p.,  p^  and  P.  and  P,  are  the  incident  and  final  electron 
and  proton  four-momenta  respectively  and  c^q  is  the  four- 
momentum  transferred  in  the  scattering  process,  i.e., 

*2c2q2=  (p±-Pf)2=  {(p-Pf)2-  (Ei-Ef)2}=  (P±-Pf)2. 

In  order  to  obtain  the  cross  section  for  this  elastic 
scattering  process  we  make  use  of    S-matrix  formulism.   The 
S-matrix  is  a  matrix  connecting  the  initial  and  final  states 
of  the  process  shown  above.   The  electromagnetic  interaction 
responsible  for  electron  scattering  may  be  treated  by  per- 
turbation theory.   The  first  approximation  of  the  S-matrix 
is  given  by  an  integral  over  the  interaction  Hamiltonian, 
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$dt  H±(t)f 


where  H. (t)  may  be  written  as  a  three  dimensional  integral 
over  a  Hamiltonian  density  H. (r,t) 

C  3   ~   -• 
H.  (t)  -  \ d  r  H.  (r,t)  . 
1       u      1 

Since  the  electrons  can  be  described  by  a  second  quantized 

Dirac  field  and  the  proton's  electromagnetic  field  may  be 

described  by  a  vector  potential,  the  interaction  Hamiltonian 

density  can  then  be  written  as 

H.(x)  =  -  j^xjA^x), 

where  x  is  the  four  dimensional  space  time  vector.   The  four- 
current  of  the  electron  can  be  written  in  terms  of  the  four 
component  state  vector  and  Dirac  matrices  as 

J®(x)  =  -  ie  Tf(x)  y^.(x),  with  Y(x)  =  Y+y4. 
The  matrix  element  of  this  scattering,  Hf . ,  is  then 

Hf.  =  -  ie  $Tf  (x)y^\.(x)A^(x)d3r  (II-l) 

giving  the  interaction  of  the  electron  current  (from  its 
initial  state  to  its  final  state)  with  the  electromagnetic 
field  of  the  proton.   By  making  use  of  the  first  Born  approxi- 
mation, the  electromagnetic  field  due  to  the  proton's  transi- 
tion current  from  the  initial  to  the  final  plane  wave  state 
can  be  obtained  by  solving  the  inhomogeneous  wave  equation 
2 


c   ot 
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where  the  expectation  value  of  the  proton's  current  operator 
is  the  transition  current  creating  the  potential  A  (x) .  The 
solution  is 

Vx)  =  *i ie  Vx,ViWi 

q 

where  9(x)  is  the  state  vector  of  the  proton.   Substituting 
this  result,  into  Equation  II-l  the  matrix  element  is 

Hfi=  -  4Tie2^Tf  (X)y  Tj  (x)  ■^?f(x)yifc«P±(x)d3x.         (II-2) 

q 

Because  of  the  use  of  initial  and  final  state  plane  waves, 
the  integration  over  space  results  in  the  delta-function 

^  (Pf+Pf  ~  P"P-)#  expressing  momentum  conservation.   It  is 

2 

interesting  to  note  that  the  1/q   term  can  be  considered  to 

be  with  either  source  current.   Thus  the  proton  may  be  con- 
sidered as  being  scattered  from  the  electron's  field  or  the 
electron  may  be  considered  to  be  scattered  from  the  field  of 
the  proton  as  it  was  above. 

The  cross  section  may  now  be  obtained  by  calculating 

4^  =  2       T~~  (H-3) 

d"   final   flux 
states 

and  averaging  over  the  initial  spin  states,  where 

"fi  =  3r  lHfil6(Ef-E.) 

and  the  flux  is  the  incoming  flux  of  electrons.   The  sum 
over  final  states  is  really  an  integral  over  all  possible 
proton  momenta  and  over  that  range  of  electron  momenta  which 
are  within  the  solid  angle  of  the  detector.   These  integrations 
over  momentum  space  and  the  averaging  over  the  initial  spins 
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and  summing  over  the  final  spin  states  are  tedious  but  then 
completed  Equation  II-3  yields 


do\  (1    e4  cos28/2\/,  ,   g2    .   2  8V 1 "\   ,,.,  A  v 

dZL.  =UlT  sTn^9/2Ai+  2n^"  tan   2  A 2ETsin^Q/2^  (II"4) 

Dirac       i       '         p  i '        ' 

proton  "    mpC2 

where  E.  is  the  energy  of  the  incident  electron  and  9  is  the 
angle  through  which  it  is  scattered.   Both  are  measured  in 
the  system  where,  the  proton  is  initially  at  rest  and 

2  ?  2        '     9 

C  ft  q   =  4E.Efsin  8/2 

Ef  is  the  final  electron  energy  given  by 


E. 


2E_. 

mc 


1  +  — |  sin  0/2 


2 

In  the  limit  that  m  c  /E  approaches  zero  the  first 

enclosed  term  on  the  right  of  Equation  II-4  was  calculated  by 
Mott  [5]  in  1929,  using  the  Dirac   equation,  to  be  the  cross 
section  of  an  electron  scattered  from  a  coulomb  potential; 
in  the  case  of  a  proton  that  physically  means  that  the  proton 
is  a  point  charge  and  that  it  has  zero  spin  and  an  infinite 
mass,  i.e.,  it  doesn't  recoil.   Let  us  call  that  the  histori- 
cal Mott  cross  section.   The  final  enclosed  term  on  the  right 
is  the  recoil  term  which  accounts  for  both  the  kinematics  of 
the  recoil  and  the  fields  associated  with  the  accelerating 
proton.   When  the  historical  Mott  cross  section  is  multiplied 
by  this  recoil  term  it  is  generally  called  the  Mott  cross 
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section.   That  nomenclature  is  followed  in  this  work.   The 

2 
tan  9/2  part  of  the  second  enclosed  term  accounts  for  the 

scattering  contributed  by  the  Dirac  magnetic  moment  of  the 

proton  while  the  1  in  that  term  accounts  for  the  point  charge 

scattering.   Equation  II-4  thus  provides  the  cross  section 

for  an  electron  scattering  from  a  point  proton  with  a  point 

Dirac  magnetic  moment.   However,  there  are  two  corrections 

which  must  be  made  to  Equation  II-4  in  order  that  the  theo- 

rectically  predicted  cross  section  of  a  proton  agrees  with 

measurements. 

The  first  thing  which  is  known  to  be  wrong  with  Equation 

II-4  is  that  the  magnetic  moment  of  the  proton  is  not  equal 

to  one  nuclear  magneton  as  is  predicted  by  the  Dirac  equation, 

but  is  measured  to  be  M-   =  +  2.7927  5  +  0.0003  nuclear  mag- 

P  - 

netons.   Also  the  prediction  of  the  Dirac  equation  of  zero 
magnetic  moment  for  the  neutron  compares  badly  with  the 
measured  neutron  magnetic  moment  of  M-   =  -1.9135  +  0.0003 
nuclear  magnetons. 

Pauli  [23]  has  shown  that  the  Dirac  equation  for  a 
charged  particle  interacting  with  an  electromagnetic  field 
can  be  modified  so  as  to  represent  a  particle  with  an  arbi- 
trary magnetic  moment  by  addding  the  operator 

V        l-L 

to  the  equation.   Foldy  [24]  has  verified  that  such  a  term 
satifies  the  restrictions  of  being  Lorentz  covariant  and  gauge 
invariant  and  has  shown  that  the  general  interaction  terms 
form  an  infinite  series  involving  arbitrarily  high  derivitives 
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of  the  electromagnetic  potential  evaluated  at  the  position 
of  the  particle.   The  series  of  coefficients  of  these  terms 
can  be  interpreted  as  a  series  of  moments  of  the  charge  and 
current  distributions  associated  with  the  particle.   We  will 
return  to  this  expansion  later  in  this  section. 

Because  the  four  potential,  h    ,■   can  be  regarded  as  the 
potential  produced  by  the  electron,  the  Pauli  term  can  be 
shown  to  reduce  to 


3C 


_e_ 


CTn.s,  (P.    "  P- 


2m  2   ^  ^v   *fv 
c 

which  is  added  to  the  y,,  term  of  the  nuclear  current.   Here 

K  =  2.7  9-  1.0  =  1.79  is  the  anomalous  magnetic  moment  of  the 
P 

proton  in  units  of  nuclear  magnetons  and 

%v  =   ~    2^"Vv~  7v7^  ' 

When  the  Pauli  term  is  inserted  into  the  matrix  element  of 
Equation  11.-2,  the  cross  section  for  a  point  Dirac  particle 
with  its  normal  point  magnetic  moment  and  now  also  its 
anomalous  point  magnetic  moment  is  obtained: 

®    •  t  =  (fS)M  J1+  4&*  fc("Va*-a<t>*  &}. 

point        Mott     **u  **  t  v 

proton 


where 


(£)..  ..-(^!^)( 


1  Mott    4E2  sin49/2/V    2Esin2Q/2 

1+      2 
m  c 
P 
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The  second  change  which  must  be  included  in  the  description 
of  scattering  from  a  real  proton  will  now  be  considered.   In 
1955  Hofstadter  and  McAllister  [22]  measured  the  cross  section 
of  the  proton  and  showed  that  the  structure  of  a  proton  is 
more  conrplicated  than  a  point  charge  and  point  magnetic  moment; 
the  proton  scattered  as  though  it  had  a  charge  and  magnetic 
moment  "cloud"  spread  out  to  an  rms  radius  of  about  0.7  fm. 
Subsequent  work  (1960-61)  on  the  proton  at  Stanford  [25,11] 
showed  that  the  magnetic  moment  associated  with  the  Pauli  term 
had  a  softer  or  more  spread  out  distribution  with  an  rms  radius 
of  about  1  fm,  while  the  distribution  of  the  charge  had  an  rms 
radius  of  about  0.8  fm.   These  rms  radii  were  calculated  on 

the  basis  of  Equation  1-3  by  assigning  proton  form  factors 

2 
F,  (q  )  to  the  Dirac  charge  and  magnetic  moment  distribution 

2 

and  F„  (q  )  to  the  Pauli  magnetic  moment  distribution.   The 

addition  of  the.  anomalous  magnetic  moment  to  the  interaction 
and  now  the  provision  permitting  the  Dirac  and  Pauli  components 
to  have  spatial  extent  changes  the  proton  contribution  in  the 
matrix  element  of  Equation  I I -2  from 


to 


>f(x)y^.<x) 


The  elastic  electron-proton  cross  section  with  this 
modified  current  matrix  element  has  been  calculated  by 
Rosenbluth  [26]  to  be 

and  is  called  the  Rosenbluth  cross  section. 
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B.   NEUTRON  FORM  FACTORS 

The  neutron  has  zero  total  charge.   Both  the  electron  and 
neutron  have  magnetic  moments,  however,  and  thus  they  can 
interact  through  the  magnetic  dipole-dipole  interaction. 
Also,  an  incident  electron  has  a  magnetic  field  associated 
with  it  because  of  its  charge  and  this  will  interact  with 
magnetic  moment  of  the  neutron.   These  interactions  are  well 
understood  and  must  be  accounted  for  in  the  calculation  of 
the  electron  neutron  scattering  formula. 

There  remains  a  third  interaction  which  is  generally  called 
the  neutron-electron  interaction  in  the  literature.   This  in- 
teraction is  expected  to  be  both  spin  and  velocity  independent 
if  there  is  a  separation  of  electrical  charge  in  the  neutron. 
Thus  an  electron,  or  any  charged  particle,  penetrating  the 
extended  charge  distribution  of  the  neutron  will  be  subjected 
to  electrostatic  forces.   These  may  arise  from  two  sources: 
1),  the  intrinsic  charge  separation  in  the  neutron  analogous 
to  the  spreading  out  of  charge  in  the  Dirac  proton  and  evi- 
denced by  the  form  factor  F^   analogous  to  F    of  the  proton, 
and  2) ,  the  anomalous  magnetic  moment  of  the  neutron,  des- 
cribed by  a  form  factor  F~  ,  analogous  to  F~   for  the  proton, 
if  the  anomalous  magnetic  moment  has  a  spatial  extent. 
Foldy  [27,28]  has  shown  that  since  the  neutron  satisfies  the 
Dirac  equation,  and  since  it  possesses  an  anomalous  magnetic 
moment,  then  some  charge  separation  is  to  be  expected.   He 
compared  the  amount  of  charge  separation  predicted  because  of 
the  neutron's  anomalous  magnetic  moment  with  that  charge 
separation  inferred  by  the  very  accurate  thermal 
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neutron-electron  measurements  of  the  neutron's  attractive 
potential  for  the  electron  made  by  several  groups  [16,17,29,30] 
The  predicted  and  inferred  charge  separations  agreed  within 
experimental  error,  implying  that  the  so-called  intrinsic  or 
Dirac  separation  is  not  present.   Since  this  requirement  for 
charge  separation  within  the  neutron  simply  because  it  has 
an  anomalous  magnetic  moment  is  not  intuitive,  this  will  be 
examined  in  more  detail  below,  as  will  the  results  of  the 
experiments  that  measured  the  neutron's  attractive  potential. 

Several  methods  of  measuring  the  attractive  potential 
have  been  used.   Only  the  first  [16]  of  these  is  outlined 
here.   If  one  measures  the  coherent  scattering  of  neutrons 
by  spinless  (in  order  to  remove  effects  due  to  magnetic  in- 
teractions) atoms  containing  a  nucleus  and  Z  electrons,  then 
the  coherent  scattering  amplitude  for  neutrons  whose  wave- 
length is  long  compared  to  the  size  of  the  atom  will  be  the 
algebraic  sum  of  the  scattering  amplitudes  due  to  the  nucleus 
and  the  Z  electrons.   By  decreasing  the  neutron  wavelength 
until  it  is  comparable  to  the  size  of  the  electronic  cloud  of 
the  atom,  destructive  interference  between  the  scattering  from 
the  individual  electrons  within  the  cloud  of  electrons  can  be 
achieved  and  the  total  coherent  scattering  cross  section  can 
be  made  to  approach  that  of  the  nucleus  alone.   By  separating 
these  scattering  effects,  one  obtains  the  so-called  neutron- 
electron  interaction.   The  result  of  the  thermal  neutron- 
electron  experiments  is  expressed  as  the  volume  integral  of 
the  potential  of  the  neutron-electron  interaction.   Plots  of 
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the  neutron  (and  the  proton  for  comparison  purposes)  charge 
density  and  the  interaction  potential  are  given  in  Figure 
II-l  [31] . 


pH         NEUTRON 


PROTON 


Figure  II-l 


Neutron  and  Proton  Charge  Densities  and  Interation 
Potentials  Between  These  Particles  and  an  Electron 


The  neutron-electron  potential  is  denoted  as  6\7  ;  that 

e  n 

is,  it  is  the  difference  between  what  the  neutron  potential 
would  be  if  it  had  no  charge  distribution  and  what  it  would 
be  if  it  had  a  positive  center  and  a  negative  outer  region. 
This  concept  of  a  positive  center-  and  a  negative  outer  region 
is  justified  by  the  results  of  the  thermal  neutron-electron 
scattering  experiments  and  is  also  consistent  with  meson 
theory.   The  volume  integral  of  this  charge  distribution 
does  not  permit  an  independent  determination  of  the  depth 

and  volume  of  the  potential.   It  does  however  permit  the  de- 

2 

termination  of  the  r     of  the  charge  distribution  by  the 

rms  3  J 

following  calculation. 
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Let 


r2    =  \r2p(r)d3r 
rms    J 


and  then  use  Poisson's  equation  to  obtain 

2  \6vd3r.  (II-6) 


rms    4TTe2 

However,  following  Fermi,  all  experiments  have  reported 
results  in  the  form  of  a  potential  V  which  was  the  average 
V  which  the  electron  sensed.   This  required  that  some  aver- 
age radius  of  this  potential  be  adopted.   By  convention,  that 
radius  is  the  classical  electron  radius  (2.81794  fm) ,  which 
implies  that  classically  the  bound  electrons  in  the  atom 
were  assumed  to  be  interacting  almost  uniformly  over  the 

neutron,  because  as  will  be  seen,  the  neutrons  largest  r 

^      rms 

is  less  than  1  fm.   That  is,  the  electrons  are  not  just 
sensing  the  outer  charge  layer  as  might  be  implied  by  the 
classical  picture  of  a  point  electron  of  low  momentum  pene- 
trating only  the  outer  shell  of  the  neutron  charge  distri- 
bution.  This  is  perhaps  best  seen  from  the  virtual  photon 
scattering  picture.   The  wavelength  of  the  virtual  photon  is 
very  long  compared  to  the  neutron  diameter  and  the  electro- 
magnetic fields  will  be  effectively  constant  at  any  time  over 
the  extent  of  the  neutron  with  the  greatest  difference  in 
these  fields  occurring  between  the  extreme  edges  of  the 
neutron.   Thus,  while  it  is  the  charge  distribution  near  the 
edges  of  the  neutron  which  show  the  deviation  of  the  neutron 
from  a  particle  with  zero  charge  density  everywhere,  it.  is  the 
average  charge  density,  and  through  Poisson's  equation  the 
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average  potential  within  the  neutron,  that  is  measured  in 
these  thermal  neutron  experiments.   Using  this  average  po- 
tential the  volume  integral  of  the  interaction  potential  is 
then  written 

r  r  XT      2 

n      e      ~>  p      e         ,  ,         ~  .         V   e         ~ 

\        Vd   r  =   V    \        4nr    dr   =  ^  m:   V  =  ~  tt  -2~-  r  (II-7) 

J  oj  3eo3  2e 

o  o  rnc 

where  r   is  the  classical  electron  radius. 
e 

Since  Equations  II-6  and  II-7  are  simply  different  ways 
of  performing  the  volume -integral  they  must  be  equal,  there- 
fore we  have: 

2      2Vo  ,  2. 

r    =  — £  (r  ) 
rms      2         e 
mc 

thus 

o      2V   /  2  N2 

r2   =  -f  (— y) 

rms      2  \   2/ 

mc    mc 

where  V   is  put  in  the  conventional  units  of  million  electron 


volts.   Then 

4  =  3i-°78  vo  (fm)2  •  (ii-8) 

The  presently  accepted  value  of  V   is  [4] : 

V   =  3720  +  90  electron  volts  .  (II-9) 

o        — 

This  value  of  V  when  substituted  into  Equation  II -8  gives 

2     '  2 

r     charge  of  neutron  =  -  (0.34  fm)  ,  the  negative  sign 

showing,  since  the  charge  on  the  outer  shell  is  weighted  by 

2 
r  more  heavily  than  the  inner  shell,  that  there  is  a  charge 

separation  in  the  neutron,  and  that  the  outer  regions  of  the 

neutron  are  negative. 
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It  was  stated  above  that  the  two  facts  that  the  neutron 
satisfies  the  Dirac  equation  and  that  it  has  an  anomalous 
magnetic  moment  are  sufficient  to  account  completely  for  the 
measured  value  of  the  potential  V  [27, 31, 32].   Let  us  now  see 
why  that  statement  is  correct. 

The  relativistic  Hamiltonian  for  a  neutron  with  the 
anomalous  magnetic  moment  as  attributed  by  Pauli  can  be 
written  as 


H  =  /3m  +  a. p-Kn (e/2m)  [£cr.H  -.ijSa-Ej 


where  K      is  the  anomalous  (and  total)  magnetic  moment  of  the 
neutron  measured  in  nuclear  magnetons.   This  can  be  reduced 
to  the  corresponding  nonrelativistic  Hamiltonian  by  the  method 
of  Foldy  and  Wouthuysen  [33]  which  yields 

H  =  £M  +  £p2/2M  -  K    (-\)/S  div  E 
11  4M^ 


Kn{-^2)P°-  [P  x  E  -  E  x  P ] 


4M 

2 

with  terms  of  order  greater  than  (1/M)  ~  being  neglected. 

For  the  coulomb  field  of  a  point  electron  located  at  x 
the  above  Hamiltonian  becomes 


H  =  £M+(|3p2/2M)+  4TTJCn(e2iS/4M2)6(x~xe)+... 

where  the  term  containing  the  delta  function  is  in  the  form 
of  the  electron-neutron  interaction.   Equating  the  volume 
integral  of  the  neutron-electron  interaction  potential  to 
this  term  one  has  [27] 

Vo  =  nKn(e2/M2)  [  (4tt/3)  (r^)]"1  =  -  3900  eV. 
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This  compares  very  well  with  the  experimental  result  given  in 

Equation  I 1-9.   Since  the  neutron  charge  distribution  can  be 

completely  accounted  for  by  its  anomalous  magnetic  moment  this 

means  that  there  is,  within  experimental  errors,  no  Dirac 

particle  intrinsic  charge  distribution  unless  relativistic 

contributions  fortuitously  cancel  this  Dirac  intrinsic  charge 

distribution. 

The  distribution  of  the  neutron  anomalous  magnetic  moment 

was  measured  by  Yearian  and  Hofstadter  and  Bumiller  [34,35] 

in  1958.   They  used  Rosenbluth's  proton  cross  section  as  given 

in  Equation  II-5  with  modifications  to  obtain  the  neutron 

cross  section.   Based  principally  on  the  thermal  neutron  work 

discussed  above,  they  set  F    =  0  to  describe  the  neutron's 

lack  of  intrinsic  charge  and  magnetic  moment  distribution  and 

allowed  F„   to  describe  its  anomalous  moment.   The  neutron 
2n 


cross  section  was  then 


32  „  (32)        F2    K2  -4[2tan2  I  +  1 

d"    WHott  2n   n  4M2L       2 


Electrons  v/ere  scattered  from  a  deuterium  gas  target  at  suffi- 
ciently high  energies  and  large  angles  to  enhance  the  magnetic 
scattering  from  the  nucleons.  By  using  the  rather  good  approx- 
imation for  inelastic  scattering  from  a  deuteron,  i.e., 

IN 
Oj   =  cr   +  or    they  were  able  to  obtain  the  r     of  the  ano- 
d     p     T}     ■*  rms 

malous  magnetic  moment  of  the  neutron  and  concluded  that  it 
was  between  0.8  and  0.9  Fermi. 
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C.   CHARGE,  MAGNETIC  AND  ISOTOPIC  FORM  FACTORS 

Thus  far  the  four  form  factors  used  to  describe  the  dis- 
tributions of  charge  and  magnetic  moment  of  the  Dirac  particle 
and  the  Pauli  moment  have  been  introduced  in  the  historical 
manner  in  which  they  evolved.   They  may  be  summarized  as 
follows : 

F,   is  the  form  factor  associated  with  the  distribution 
lp 

of  charge  of  the  Dirac  proton  and  also  with  the  distribution 

of  the  intrinsic  (Dirac)  magnetic  moment.   In  the  static 

(q2  =  0)  limit  F   (0)  =  1. 

F~   is  the  form  factor  associated  with  the  anomalous 
2p 

(Pauli)  magnetic  moment  distribution  of  the  proton.   In  the 

static  limit  F0  (0)  =  K  =  1.79. 
2pv      p 

F,   is  the  form  factor  associated  with  the  distribution  of 
In 

charge  of  the  Dirac  neutron  and  also  with  the  distribution  of 

the  intrinsic  (Dirac)  magnetic  moment  of  the  neutron.   In  the 

static  limit  Fn  (0)  =  0. 
In  v  ' 

F„   is  the  form  factor  associated  with  the  anomalous 
(Pauli)  magnetic  moment  distribution  of  the  neutron.   In  the 
static  limit  F2n(0)  =   \   =    1.91. 

It  has  been  shown  above  that  the  anomalous  (Pauli)  mag- 
netic moment  of  the  neutron  causes  some  charge  distribution 
within  the  neutron.   Therefore  the  form  factor  F~   not  only 
tells  what  the  anomalous  magnetic  moment  distribution  is,  but 
tells  something  about  the  charge  distribution. 

The  work  of  Yennie,  et  al.  (1957  [31],  Ernst  et  al .  (1960) 
[8]  and  Hand  et  al.  (1963) [19]  introduced  linear  combinations 
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of  the  four  Dirac  and  Pauli  form  factors  for  both  theoretical 

and  experimental  reasons.   These  combinations  were  the 

following: 

2 

G„  =  F, %- r  ^  Fn 

E     lp   „  2  2   p  2p 

p     *   4m  c       * 

G    =  Ft   +  K  F 
mp     lp     p  2p 

G    =  F    — g —  k  f 

^E„   *ln   ,22   n^n 
n       4m  c 

G   =  F,  +  ^  F   . 
m     In   n  2n 

In  Ref.  19  it  is  pointed  out  that  these  linear  combinations 
correspond  to  zero,  G  ,  and  one,  G  ,  unit  of  angular  momentum 
transferred  along  the  direction  of  the  virtual  photon  ex- 
changed in  the  scattering  process.   Reference  [8]  added  more 
physical  content  to  these  linear  combinations  of  F,  and  F~  by 
showing  that  the  rms  radii  of  and  charge  and  magnetic  moment 
spatial  distribution  are  derived  from  Equation  1-3  by 


[6AGe]2       and  t6 


d_ 
2   m 


dq  dq 

2 

at  q   =0.   For  this  reason  the  descriptive  names  for  Gp  and 

G   are  the  "charge"  and  "magnetic"  form  factors.   Additionally, 
Ref.  19  points  out  that  the  cross  terms  of  the  charge  and  mag- 
netic form  factors  do  not  interfere  with  each  other  in  the 
Rosenbluth  cross  section  when  it  is  written  in  terms  of  G„ 

and  Gx„,  while  the  Dirac  and  Pauli  terms  are  mixed  into  the 

M 

Rosenbluth  cross  section  as  products  of  each  other,  as  can 
been  seen  by  examining  Formula  II-5.   Written  in  terms  of  Gp 
and  G  ,  the  Rosenbluth  cross  section  for  the  proton  becomes 
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Mott  1+  3__   4mp   1+  a__  P 

4m  4m 

P  P 


with  an  analogous  equation  for  the  neutron  cross  section. 

G   and  G   are  obtained  experimentally  by  measuring  the  cross 

section  at  at  least  two  different  scattering  angles,  at  the  same 

2  2 

q  .   From  the  intercept  of  a  plot  of  G   versus  the  coefficient 

2 

of  G„„   one  obtains  G„  and  from  the  slope  one  obtains  G.,  . 
ME  c  M 

P        2        P  P 

The  static  (q  =0)  values  of  these  charge  and  magnetic  form 

factors  are: 

G   (q2  =  0)  =  1 

EP 

GM  (q2  =  0)  =  2'79  (II-ll) 

P 

GE  (q2  =  0)  =  0 

n 

GM  (q2  =  0)  =  -  1.91. 
n 

We  now  consider  what  the  slope  of  G    tells  about  the 

n 
charge  distribution  of  the  neutron  by  studying  the  primary 

scattering  event  of  interest  in  this  work,  the  e-d  scattering.. 

This  process  may  be  described  by  the  Feynman  diagram 
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where  e.  and  ef  are  the  incident  and  scattered  electron 
states,  d.  and  df  are  the  target  and  recoiling  deuteron 
states  and  q  is  the  four -momentum  transfer.   The  deuteron  has 
a  diameter  of  about  4  fm;  in  this  work,  with  electron  energies 
of  about  35  to  100  MeV,  the  virtual  photons  had  wavelengths 
in  the  range  of  10  to  30  fm.   With  these  long  wavelengths 
compared  to  the  extent  of  the  deuteron  and  neutron,  only  the 
outer  charge  distributions  are  measured.   It  is  important  to 
note  that  the  thermal  neutron  experiments  [18,4,16,17]  give  a 

positive  slope  for  — -j  G    as  do  two  of  the  e-d  experiments 

dq  12 

[20,21].   This  means  that  the  quantity  -  -y  r  is  positive, 

i.e.,  that  the  outer  part  of  the  neutron  is  negative.   This 

is  so  because  the  square  of  the  rms  radius  is  obtained  from 

the  expression 

CO 

r     =  \   r  (4nr  p(r)dr) 
rms    J     v     t-  \    i       i 
o 

CO 

where    the    total    charge   of   the   neutron   =    \      4vrr    p(r)dr  =    0; 

o 

shows  that  the  chare 
rms 

of  the  neutron  must  be  more  negative  at  larger  radii. 


2 
thus  a  negative  r     shows  that  the  charge  density,  p(r), 


Historically,  theoretical  attempts  have  been  made  to  explain 
nucleon  charge  distributions  in  terms  of  a  third  set  of  nucleon 
form  factors  which  are  called  the  isotopic  form  factors.   These 
isotopic  form  factors  are  useful  because  the  structure  of  nu- 
cleons  is  determined  by  the  strong  interaction  and  calculations 
can  be  arranged  in  such  a  way  that  one  can  take  advantage  of 
the  conservation  of  isospin.   As  the  form  of  the  nuclear  elec- 
tric current  operator  [36] 

1+  T 
J^(x)  =  e«p.(x)y^  — 2"-^-  VM 
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shows,  the  current  has  mixed  transformation  properties  under 
rotations  in  isospin  space;  part  of  it  transforms  as  a  scalar 
and  part  as  a  third  component  of  an  isovector.   As  a  conse- 
quence of  these  mixed  transformation  properties,  the  proton 
and  neutron  current  elements  may  be  expressed  in  terms  of  an 
isovector  and  isoscalar  current,  i.e., 

j£(x)  =  jS(x)  -  jV(x). 

The  electric  and  magnetic  isotopic  form  factors  then 
enter_  the  cross  section  calculation  as  the  coefficients  of 
the  scalar  and  vector  current  operators  in  the  matrix  element. 
The  isotopic  form  factors,  where  the  superscript  denotes 
scalar  or  vector,  are: 

4  -  *  h  +  ge  ]  7^T h 

p     n   - 


ge  =  *  k  -  s  ]  7^ h 

p    n   n 

p    n   ^  ^ 

G^   =    ^    [gm    -  GM    1      T~TZ*k     (M<    -H    ) 

M         ^Lm  MJq^O       ^     v    p      n 


II-12) 


There  have  been  attempts  to  give  the  isotopic  form  factors 
further  meaning.   While  these  attempts  have  not  been  completely 
successful,  it  is  useful  to  consider  them  briefly.   A  state  of 
n  tt  mesons,  produced  by  a  photon,  must  have  zero  total  charge 
and  be  odd  under  charge  conjugation.   Based  on  this,  Drell  and 
Zachariasen  [37]  have  shown  that  only  an  even  number  of  pions 


37 


contribute  to  the  isotopic  vector  form  factor  and  that  only 
an  odd  number  of  x^ions  can  contribute  to  the  isotopic  scalar 
form  factors.   They  also  have  shown  that  charge  conjugation 
does  not  permit  a  single  pion  to  be  produced  since  the  photon 
is  odd  under  charge  conjugation  and  the  n  ,  the  only  ti  which 
conserves  charge,  is  even  under  charge  conjugation.   The 
possible  intermediate  states  then  include: 


Figure  II-2 
Possible  Intermediate  Photon-Proton  States 

In  the  diagram  in  Figure  II-2  the  possible  existence  of 
pion-pion  forces  were  disregarded;  but  we  know  that  there  are 
resonances  in  the  two-pion  and  three-pion  systems,  and  we 
expect  that  such  resonances  may  have  an  effect  on  the  elec- 
tromagnetic structure  of  the  nucleons.  Chew  et.  al.  [38]  and 
Nambu  [39]  have  shown  that  such  propagators  as  depicted  in 
Figure  II-2  (b)  and  (c)  do  in  fact  contribute  to  the  form 
factors  in  such  a  way  that  they  have  poles  at  the  masses  of 
the  p  ,    co  and  cp  .   The  rho  (isospin  =  1,  spin  =  1)  contributes 
to  the  isovector  and  the  omega  and  phi  (isospin  =  0,  spin  =  1) 
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contribute  to  the  isoscalar  form  factors.   On  the  assumption 
that  virtual  clouds  of  mutually  strongly  interacting  bound 
pions  are  responsible  for  the  electromagnetic  properties, 
these  clouds  show  up  as  resonances  in  dispersion  theoretical 
expressions  for  the  form  factors,  for  instance 

2S  .  .2 


q2)  =  M   ^^^F-  (H-13) 

i-i-i       rr   '    _r-f 


m    q-  -q 


where  m  corresponds  to  the  threshold  for  resonance  creation. 

Furthermore,  if  one  assumes  that  these  vector  mesons  (spin  =1) 

2 

appear  as  sharp  peaks  in  the  spectra  of  weight  functions  g(q'  ), 

i.e.,  neglecting  the  widths  (P(100  MeV)  ,  to(10  MeV)  ,  9(3  MeV)  , 
Equation  11-13  becomes 


G|(q2)  =|[(l-Sel)+— -Si],  (11-14, 

i+  a 

which  is  the  form  factor  expression  which  Clementel  and 
Villi  [40]  have  proposed,  where  S  ,  is  a  constant.   The  term 
1-S  -.  represents  the  possibility  of  high  mass  resonances, 
i.e.,  it  represents  those  resonances  which  have  such  high 
masses  that  their  poles  may  be  replaced  by  constants.   M   is 
the  mass  of  the  scalar  resonances.   There  is  one  term  for 
each  resonance.   Similar  expressions  for  G  ,  G   and  G  were 
likewise  developed  and  were  used  by  deVries  [10]  to  fit  proton 
and  neutron  data.   This  fit  is  discussed  in  Appendix  A  and 
provided  the  normalization  for  the  proton  data  obtained  in 
this  work. 
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Earlier  in  Section  II  it  was  observed  that  the  outer 
regions  of  the  neutron  have  been  measured  to  be  negative  in 
the  thermal  neutron  work.   There  are  several  theoretical 
reasons  to  believe  that  the  outer  regions  should  be  negative 
[32,38],   Among  these  is  the  fact  that  the  anomalous  magnetic 
moment  is  almost  entirely  a  vector  in  isotopic  spin  space, 
i.e.,  the  anomalous  moments  of  the  neutron  and  proton  are 
nearly  equal  in  magnitude  with  opposite  signs.   /additionally 
since  the  isovector  meson,  the  rho,  is  the  least  massive  of 
the  rho,  omega,  phi  mesons,  the  uncertainty  principle  as  well 
as  the  Yukawa  potential  indicate  that  the  outer  regions  will 
be  dominated  by  the  rho  meson. 


D.   METHOD  OF  EXTRACTING  THE  PROTON  ELECTRIC  FORM  FACTOR  FROM 
THE  PROTON  CROSS  SECTION 

The  Rosenbluth  cross  section  for  the  proton  (Eq.  11-10) 

can  be  written  as: 

proton     Mott  p      p 

2    2 

where  T  =  q  '/4M  .   This  can  be  simplified  by  using  the  scaling 

law  proposed  by  Lehmann  et.  al.  [41]  and  others  [19]  ,  which 

states  that  G 

M 
G  =  ^  .  (II-15a) 

P    P 

Schumacher  has  shown  that  we  should  not  expect  isoscalar 

form  factors  to  equal  isovector  form  factors  [42] ,  and  that  a 

set  of  scaling  laws  which  relate  the  form  factors  G    and  G 

P       n 
to  GE   and  G    are 

P        n        GM   -  Ve  +  ^nGE 
P        P       n 
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and 


M   =   n  E     p  E 

n        p    *      n 


For  our  experiments,  magnetic  scattering  generally  accounted 
for  only  a  percent  or  two  of  the  total  scattering  in  most  runs 
and  at  most  10%  of  the  scattering  in  a  few  runs;  thus  the  dif- 
ference in  these  two  scaling  laws,  which  was  at  most  a  1.0% 

2 

difference  at  0.35q  ,  then  made  a  difference  in  the  cross 

section  of  less  than  a  nominal  0.01%  up  to  at  most  a  0.05%. 
Thus  the  conventional  scaling  law  was  used  to  reduce  the  cal- 
culations.  When  the  magnetic  contribution  is  accounted  for 
by  using  the  scaling  relationship,  Equation  II-15a,  Equation 
11-15  reduces  to 

4 

h~  =  °£(e)  TkF  L1+  r^+d+r)  (2r)M.2tan2  §].  (ii-ie) 

where 

Mott 
Because  the  last  two  terms  in  the  bracket  are  due  to  the 
magnetic  form  factors  it  is  convenient  to  group  them  together 
by  defining 

Cp  =  THp+  2T^(1+T)tan2  §  .  (11-17) 

Because  the  proton  cross  section  can  be  written  as 

dcr^     r>     2 

h~  -  CTM(e»Gp-  ("-18' 

where  G   is  the  proton  form  factor,  we  can  obtain  the  proton 
electric  form  factor  by  combining  11-16,  11-17,  and  11-18  as 


(11-19) 


da? 

G2    = 

g2,-(i+t; 

=    K?G2 

=  KP     d* 

E 
P 

p<i+cp) 

M   p 

M   «S(B) 
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where 

P  _  u±n„ 

KM  ~  (1+C  ) 

is  defined  as  the  proton  magnetic  correction  term.   Thus  by 
measuring  the  proton  cross  section  and  dividing  it  by  a^(0) 
and  then  multiplying  this  quotient  by  kE  one  extracts  the 
square  of  the  proton  electric  form  factor. 


E.   THE  DEUTERON  GROUND  STATE 

Before  the  elastic  electron-deuteron  calculation  is  studied 
it  is  useful  to  consider  the  ground  state  of  the  deuteron. 
This  is  particularly  true  in  this  case  of  elastic  scattering 
because  the  deuteron  is  in  its  ground  state  both  before  and 
after  the  scattering  event.   It  is  an  experimental  fact  that 
the  deuteron  has  but  one  bound  state  and  that  is  the  ground 
state.   Since  the  discovery  of  the  deuteron  in  1932  many  of 
its  characteristics  have  been  measured  and  provide  useful 
guidance  in  interpreting  the  electron  deuteron  scattering  data. 
Some  of  the  more  important  are  discussed  below. 

1)  The  neutron  and  proton  are  bound  together  by  an  energy 
of  2.2256  +  0.0015  MeV. 

2)  The  angular  momentum  is  1ft,  and  the  parity  is  even. 

3)  The  magnetic  dipole  moment  is  n,  =  0.857  393+0.00003 
nuclear  magnetons. 

4)  The  electric  quadrupole  moment  is  0.00282  +  0.00002. 

5)  The  root-mean-square  electromagnetic  radius  of  the 
deuteron  is  1.96  +  0.02F. 

6)  The  deuteron  potential  has  a  repulsive  core  with  an  rms 
radius  of  about  0.4F. 
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7)  The  mass  of  the  dcuteron  is  187  5.58  5  MeV. 

Elementary  calculations  solve  for  the  deuteron  wave 

function  in  a  square  well  neutron  proton  potential  using  the 

measured  binding  energy.   The  solution  shows  that  the  neutron 

and  proton  spend  most  of  the  time  outside  of  the  square  well 

potential.   The  ground  state  is  mostly  a  IS  state,  i.e.,  an 

1=0    state,  which  is  indicated  not  only  from  general  energy 

considerations  [43] ,  but  from  the  fact  that  the  algebraic  sum 

of  the  proton  and  neutron  magnetic  moments  nearly  equals  the 

deuteron  magnetic  moment.   This  indicates  that  the  intrinsic 

spins  of  these  particles  are  aligned  and  that  there  is  no 

orbital  angular  momentum.   If  the  nuclear  force  were  spin 

independent,  one  should  expect  to  observe  a  singlet  state, 

1    with  the  same  energy.   However,  no  bound  state  with  J  =  0 

0 
has  been  found  in  the  deuteron.   This  indicates  that  the 

neutron-proton  force  is  stronger  when  the  spins  are  parallel 

than  when  they  are  antiparallel ;  however,  this  does  not  mean 

that  the  antiparallel  or  singlet  force  is  zero.   In  fact,  the 

singlet  force  is  almost  strong  enough  to  produce  a  bound  state. 

Although  the  above  description  of  the  deuteron  is  accurate 

to  first  approximation,  a  better  description  is  possible. 

First,  one  notes  that  the  algebraic  sum  of  the  proton  and 

neutron  magnetic  moment  is  0.87925  +  0.0004  while  the  measured 

magnetic  moment  of  the  deuteron  is  0.857393  +  0.00003  nm. 

This  difference,  about  2^f/o,    is  larger  than  the  experimental 

error.   The  second  problem  with  the  description  given  above  is 

that  a  pure  S  state  implies  that  the  quadrupole  moment  is  zero 

since  the  wave  function  is  spherically  symmetric,  while  the 
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deuteron  in  fact  has  a  positive  quadrupole  moment.   The  only 
other  state  in  which  two  particles  give  the  same  angular  mo- 
mentum and  parity  is  the  ^D,  state  (&   =  2,  S  =  1,  J  =  1) . 
The  conclusion  is  that  the  total  wave  function  of  the  deuteron 
must  have  the  form  Y  -  a  Y   +  a  Y   where  the  subscripts  indi- 
cate S  state  and  D  state  and  the  coefficients  satisfy 

2     2 

aq  +  a   =  1   in  order  to  normalize  the  total  wave  function. 

The  exact  value  of  a   that  produces  a  quadrupole  moment  in 
agreement  with  the  experimental  value  depends  on  the  radial 
part  of  the  function    Y  ,  which  in  turn  depends  upon  the  de- 
tails of  the  radial  dependence  of  the  neutron-proton  potential. 

2 
Values  of  the  D  state  probability,  P   =  a  ,  in  the  range  of 

3  to  7%  have  been  obtained  from  calculations  using  different 

2 

potentials  [44] .   A  mixed  deuteron  ground  state  with  a   w  4% 

produces  an  electric  quadrupole  moment  and  a  magnetic  dipole 
moment,  both  in  qualitative  agreement  with  observations.   It 
is  one  of  the  objects  of  this  work  to  help  determine  which 
wave  function  best  fits  these  values  and  from  this  determine 
more  properties  of  the  nuclear  potentieil.   Different  wave 
functions  which  are  considered  in  this  work  have  different 
percentages  of  D  state.   However  it  is  possible  to  keep  the 
same  percentage  of  D  state  and  to  change  the  quadrupole  moment 
by  changing  the  range  of  the  tensor  potential.   This  potential 
cannot  be  spherical  as  that  produces  a  pure  S  state.   The 
positive  quadrupole  moment  indicates  that  the  two  nucleons 
prefer  to  line  up  with  their  spin  vectors  along  the  same  line 
rather  than  side  by  side.   The  tensor  potential  which  is  a 
function  of  the  angle  between  the  spin  vector  of  the  neutron 
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and  the  proton  and  the  radius  vector  r  separating  them  accounts 
for  this  preferential  alignment  of  spins  and  is  a  potential 
which  must  be  added  to  the  dominant  spherical  potential. 

Several  textbooks  [43,45,40]  cover  in  detail  all  the  con- 
cepts that  were  briefly  discussed  above.   Having  summarized 
the  deuteron  ground  state  we  now  consider  in  detail  the 
elastic  electron  deuteron  scattering  process. 

F.   ELASTIC  ELECTRON-DEUTERON  SCATTERING 

Jankus  [13]  in  1956  derived  the  cross  section  for  elastic 
electron-deuteron  scattering  in  the  first  Born  approximation. 
In  his  calculation  the  nucleons  were  represented  by  point 
charges  and  point  magnetic  moments.   The  electron  was  con- 
sidered extremely  relativisitic,  but  the  deuteron  was  treated 
nonrelativistically. 

While  the  nucleons  have  spin  one  half  and  two  form  factors, 
the  deuteron,  with  spin  1,  has  three  form  factors.   This  is 
consistent  with  the  work  of  Glaser  and  Jaksic  [47]  and  others 
[48]  which  predicts  that  the  cross  section  for  scattering  a 
relativistic  electron  from  a  potential  of  spin  J  will  contain 
2 J  +  1  form  factors. 

Jankus' s  result,  written  in  the  more  current  terms  [49], 
and  including  form  factors  for  the  proton  and  neutron,  i.e., 
permitting  Jankus 's  point  nucleons  to  have  extent,  is 

d 


d 


7T   =  aM(9){Go(q[2}+  ^G\^  +   f^(q2)[l+2  (1+77)  tan2  §]},  (11-20) 


where  17=  q  /4Md  and  GQ  (q2)  ,  G-^q2)  and  G2(q2)  are  the  charge, 
quadrupole  and  magnetic  moment  form  factors  of  the  deuteron. 
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These  are  given  by 


where 


G0  =  (GE+  GE  )CE'  (1I~21) 

P    n 


CO 

CE  =  I      [u2(r)+  w2(r)]jo(q  |)dr  (11-22) 


o 
for  the  charge  from  factor; 


G2  =  (GE  +  GE  )CQ/  (11-23) 


where 

2 


CQ  "  4.  S   [u(Dw(r)  -^  ]j2(q  |)dr        (II-24) 


Gl  =  2  J4  {(GM  +  GM  )CS+(GE  +  GE  )CJ  (lI-25; 


'Q    \[2t?  JQ  L  2f-2" 

for  the  quadrupole  form  factor,  and 

M 

1 

p      p     n        p     n 

for   the  magnetic    form   factor,    where 

C      i  T   2  #    .      w(r)~l.     ,    rN[w(r)u(r)         w    ( r )   1 .     ,      r\  ,       / ,.  -r    o ,-  \ 
CS=    i       lLu     ^)"  ~-^p0i^)+l-t±±±  +  — 2      J32(«    2)dr     (II~26) 
o  N     2 

and 

CL  =   fS      w2(r)[jo(qf)+    j2(qf)]dr.  (11-27) 


(G„,  +  G,,  )  C„  determine  the  contribution  of  the  intrinsic 
M    MS 
P     n 
magnetic  moments  of  the  proton  and  the  neutron  to  the  cross 

section  and  (G„  +  G„  )CT  determines  the  magnetic  contribution 

£i      h  JLi 

p     n 
of  the  convection  of  charge  in  the  deuteron  to  the  cross 

section. 

Here  M   and  M,  are  the  mass  of  the  proton  and  deuteron 

respectively  and  u (r)  is  the  S-state  wave  function  and  w(r) 
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is  the  D-state  wave  function  normalized  so  that 

^   [u2(r)-f  w2(r)]dr  =  1  (11-28) 

o 

and  j0(q  §0  and  j2(q  ^-)  are  spherical  Bessel  functions. 

The  static  limits  of  the  functions  above  are 

GQ(0)  =1,      G2(0)  =  M^Q, 


p 

;[ge  (0)+  Ge  (0)]  =  GE  (0)=  h,h(GM    (0)+GM  (0))=GM  (0)=%(l*  4-H^) 


CE(0)=  1,  CQ(0)=  M^Q,  Cg(0)=  1  -  |  PD  (11-29) 

and  - 

cL(0)  =  f  pd, 

Q  is  the  quadrupole  moment  of  the  deuteron,  M-   and  ^   are 

"*     *•  p      n 

the  proton  and  neutron  magnetic  moments,  and 

CO 

PD=  ^  w2(r)dr 

is  the  D-state  probability  of  the  deuteron. 

Introducing  the  nucleon  form  factors  to  Jankus's  results 
of  course  assumes  that  the  nucleons  are  not  distorted  in  the 
deuteron.   The  fact  that  these  nucleons  are  so  weakly  bound 
that  they  spend  a  sizable  portion  of  their  time  far  apart 
makes  this  assumption  of  non-distortion  plausible. 

Gross  [50]  and  Casper  [20]  have  calculated  relativistic 
corrections  to  the  electron  deuteron  scattering  problem. 
While  these  corrections  could  come  primarily  from  relativistic 
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modifications  of  (a)  the  deuteron  wave  functions  and  (b)  the 
nucleon  current,  the  latter  have  been  shown  to  be  small  by 
Gourdin  [71]  and  also  by  Gross.   The  relativistic  modifications 
of  the  deuteron  v/ave  function  depend  not  just  on  the  Lorentz- 
contraction  factor,  which  is  actually  a  contraction  in  two 
directions  [31],  but  also  on  the  fact  that  the  equal-time  wave 
function  in  the  rest  frame  is  not  the  same  as  the  equal-time 
wave  function  in  the  moving  frame  [20] .   Gross  goes  on  to 
solve  for  the  corrections  to  the  deuteron  wave  function  and 
then  obtains  workable  expressions  for  the  corrections  in  terms 
of  u(r)  and  w(r).   Gross's  expression  for  G„,  G, ,  and  G~    are 
then  much  more  complicated  than  those  of  Equation  11-21,  11-23, 
and  11-25.   For  example,  Gross's  complete  result  for  the  deu- 
teron charge  form  factor  is 

Cp  =(l-  q2/32M2)\'   (u2+w2)jn(q  §■)  dr-  -9-j  \       (u2-uru' 
E  d   o  °    Z  24M   Jo 

-  uku  +  7w  -  w  rw*  -  whw  )|_Po^q  2~^  +  ^2^q  2^  Jdr 

2   r  w 

+   q  2  \   (uru'+  wru>  +  wrw*   )  j  (q  |-)dr,  (11-30) 

16M  -.   o 
d 

where  prime  refers  to  differentiation  with  respect  to  r  and 
2 


k  =   r2 


I: — 5-  -   Of  L       a  =    (Me,    and    e   =    deuteron  binding   energy, 


"dr' 

The  first  integral  is  the  usual  Jankus  nonrelativistic  result. 
Schumacher  and  Bethe  [44]  and  Elias  et  al.  [49]  note  that 

the  relativistic  corrections  of  Casper  and  Gross  [20]  to  the 

2    2 

deuteron  form  factors,  to  order  q  /2M  ,  enter  simply  as  a 

2 
factor  (1-  T/2)   on  the  right  hand  side  of  Equation  11-20, 
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2         2 
where    T  =   q"/4M    •      Equation    11-20    can  be    further   simplified 

2 
for    small   q '   by  using    the    scaling    law    [42] 

M 

Gl   =    ^dG0   M~  (II"31) 

P 

which  is  not  a  law  but  an  approximate  fit  to  the  data  which 

can  be  derived  theoretically  from  the  nonrelativistic  quark 

model  and  SU6.   (SU6  also  predicts  M-  /M-   =1.5  which  is  not 
v         *  p'  n 

quite  true.) 

By  making  the  relativistic  correction  as  noted  by 
Schumacher  above  and  using  the  scaling  law,  Equation  11-20 
becomes 


ajJ(9){G2)(q2)  +  §U2G^(q2)+  |^G2T[l+2(l+7))tan2  §]}(!-  f) 


tisx.    - 

(11-32) 

Rewriting  this  using  Equation  11-21  and  dropping  the  quad- 

-4 
rupole  term  as  it  is  less  than  5  x  10    at  the  highest  value 

2 

of   q      reached   in    this    experiment,    we    then  have 

f~-  =    aM(9)  (GE    +GE)2ce{1+   f^r  (1+2  (1+7))  tan2    § )  }  (1-   |)2.     (11-33) 
p        n 

If  we   define   Cd  =   -|t>2  (1+2  (1+tj   tan2   -|)  )  (11-34) 

and  note    that 

i^=  ffM(9>Gd-  (II-35) 


then      2  2  2 

Gd  =  (GE  +  GE  )    CE(1+Cd} (1~  2}  '  (II-36) 

p     n 

Additionally,  if  we  define  K   as  the  deuteron  magnetic  and 


49 


rclativistic  correction  we  have: 

<r=    <Tfc7><7~T  )2-  <"-") 

a   1-  -j 

We  obtain  the  so-called  deuteron  electric  form  factor  G„ 

Ed 
from  Equation  11-36,  i.e., 

G2  d°d 

G2  .    (GE  +GE    )2C2  =  d     T2        -   GXr=Kmr  "f"      '     (II-38: 

Ed    EP  En   E  (i+cd)(i-~)2     Q  mr   mr  (£(0) 

With  both  the  deuteron  and  proton  electric  form  factors  avail- 
able we  now  obtain  the  neutron  electric  form  factor  by  first 

forming  the  ratio  of  G„   to  G   ,  i.e., 

d      £'p 


V<&       (ge  taa>  +GE   ^>ce 

d  2      = E- P ,  (H-39) 

GE  V  GE  «# 

P   F  P   F 

2       2 

where  the  q,  and  q   indicate  that  the  proton  and  neutron  in 

the  deuteron  have  received  a  four-momentum  transfer  different 
from  that  of  the  proton  which  was  the  nucleus  of  the  hydrogen 
target.   This  was  so  in  this  experiment  because  the  incident 
electron  energy  was  such  that  the  recoil  electron  from  both 
hydrogen  and  deuterium  targets  had  the  same  final  energy. 
Thus  the  transmission  efficiency  of  the  spectrometer  and  the 
detector  efficiencies  were  held  constant  for  both  targets,  but 
the  four-momentum  transfer  was  slightly  different.   Since  the 

proton  electric  form  factor  is  available  from  deVries  [10] b1 

2 
fit  we  may  scale  the  G  (a)    measured  from  hydrogen  to  the 

2  P 

GE  ^qd^  measured  in  deuterium.   Thus  the  electric  form  factor 
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of  the  proton  in  the  hydrogen  that  would  be  predicted  from 
the  hydrogen  measurement,  if  the  hydrogen  had  experienced  the 
same  four-momentum  transfer  as  did  the  proton  in  the  deuterium 
is  given  by 


P  p 


GE  (3d> 


P  * 


(11-40) 


deVrie: 


Equation  11-39  may  be  rewritten  in  terms  of  the  measured 
quantities,  the  deuteron  and  proton  cross  sections,  by  using 
Equations  11-19  and  11-38.   This  gives 


GE  (qd)+  GE  (<5> 
L  P      n - 


GE  <<# 

P   F 


(11-41) 


Now  forming  this  ratio,  i.e.,  Equation  11-41,  with  the 
same  four-momentum  transfer  for  the  proton  in  both  the 
deuterium  and  hydrogen  we  have: 


■   (q^+G*   {qh 


E  ^d' 


^   aP(0)  Kd 

d-a   m      mr 


GE  <<& 
P 


which  gives 


%£-   0d(9)  KP  /   GE 
d*.        m v     m  /  \ 


(11-42) 


deVries 


2  A 


-  1, 


(11-43) 


deVries 


where  R  is  the  experimental  ratio  of  the  deuteron  cross  section 
to  the  proton  cross  section. 
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Note  that  cp(0)  and  ¥p   wore  calculated  at  the  incident 


energy  of  the  electron  on  the  hydrogen,  which  was  higher 

than  the  energy  of  the  electron  on  the  deuteron,  at  which 

c  (9)  and  K   were  calculated, 
m   '       mr 

2 
The  neutron  electric  form  factor,  G   (q,) ,  then  is 


given  by 


v# 


oP(G)    Kd 

m  mr 

ad(9)   KP 


\l<$ 

v^ 


1.0 


deVries 


\«& 


(11-44: 


which  may  be  rewritten  a 

fa 


ge  <<#  =  1 

n 


(II-44b) 


2 

G   (q-.)  is  the  G„   which  includes  Gross's  relativistic 
E    d  E 

^n  n 

corrections.   It  appears  at  first  to  be  quite  different  from 
those  of  Stewart  and  Mader  because  they  used  the  nonrelati- 

vistic  deuteron  wave  function  to  calculate  G„   and  then  added 

n 
the  change  in  G   which  is  due  to  this  relativistic  correction 

n 
in  the  deuteron  wave  function  to  the  nonrelativistically  cal- 
culated G_  .   They  used  Gross's  equation  10  [20]  in  a  form 
E 

n  2 

which  is  practically  model  independent  within  the  range  of  q 

considered  here  [2],  i.e., 

(GE  >    -  GE  +  AG   =  GE  +  ^™p  • 
n  rel    n    *>n  n       e 

This  can  be  shown  to  be  very  close  to  the  correction  of 

Schumacher  used  in  this  work  by  considering  the  following. 

Schumacher's  relativistic  correction  enters  as  a  factor 

2 
(1  -  T/2)   on  the  right  hand  side  of  Equation  11-20.   This 
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correction  is  carried  forward  to  Equation  11-43.   If  we 
form  the  nonrelativistic  ratio  of  GE   to  GE   we  have 


Rel 


(7p(9)  Kd  r: 
_JL     m 
"d 


V# 


cra(9)  KP  /  /  G„  (qi) 
m     m  /    E   ^d 


-  1, 


(11-45) 


Define 


m 


(1-  T/2) 


If  K   were  inserted  in  Equation  11-44  in  place  of  K  we 
mr  m 

would  have 


Rel 


Since 


0-P(0)  Kdf   (GE  V 


Od(0)  R 


and 


AGE  =  I  R 
n 


p/Rel 


Qm(Q)  Km 
G>>  Km 


GE  (^d> 


CE(1-  T/2) 


-  1 


p  /non 
Rel 


UH  «4> 


%   (^> 


^  p 


'GE  V 


P 


GE  <*d> 

E 


8m 


(11-46; 


(11-47 


(11-48) 


2    2 

where  T  =  q  /4m  has  been  written  explicitly.   Since  the 

2    2 

coefficient  of  q  /8m   in  Equation  11-48  is  experimentally 

within  2  percent  of  unity  throughout  this  experiment,  then  we 

see  that  Schumacher's  method  and  Stewart's  and  Mader ' s  method 

are  in  good  agreement. 
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G.   THE  ROSENBLUTH  PLOT 

Rosenbluth  derived  the  cross  section  for  electron-proton 
scattering  based  on  the  one  photon  exchange.   Jankus  made  the 
same  assumption  in  calculating  the  deuteron  cross  section. 
Examining  Jankus' s  result,  given  in  Equation  I.I-20  ,  and,  re- 
produced below, 

fr  =  (zr)       k  <*3>+  h24  te2)+  M  ^  ^  ^ ta"2  fi  i 


/Mott 


which  can  be  written  as 


(II-48a) 


ad .  fad)       G2  .  .    (II_48b) 

it  can  be  seen  that  if  the  deuteron  cross  section  is  measured 

2  2 

at  different  angles  but  at  a  constant  q  ,  then  a  plot  of  G, 

2  G 
versus  (1+2  (1+??)  tan  ~x)    should  give  a  straight  line  with 

2    o  2  2  2   2 

G  +  q-t?  Q>      as  the  intercept  and  -o"??G,  as  the  slope.   Thus  a 

linear  plot  of  such  data  would  indicate  that  the  one  photon 

process  accurately  describes  the  phenomena.   Such  a  plot  has 

been  made  in  Section  VB. 

Equation  II~48b  cannot  be  used  directly  because  that 

equation  is  presumed  to  give  an  absolute  cross  section.   The 

absolute  deuteron  cross  section  is  obtained  through  knowledge 

of  the  hydrogen  cross  section.   Since  deVries  b'  fit  is  the 

best  fit  to  the  proton  absolute  cross  section,  this  is  used  to 

normalize  the  data.   That  is 


*&\  s  (a°P] 


do9 
dft    ^Absolute   Vd-a  ^deVries     Vd-a  /  LINAC 
where  N  is  the  normalization  factor  for  a  given  point. 
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Therefore 

\    /deVries  /  \ 


(11-50 
NPS  LINAC 


Because  the  deuteron  and  proton  cross  sections  were  measured 
under  the  same  experimental  conditions,  the  absolute  deuteron 
cross  section  is 

dsL    /Absolute      Vd^  /  NPS  LINAC 


It  follows  that 


jilted 

leVries  /  \dA    J     V^^/linac 


In  terms  of  the  experimental  ratio  R,  where 
ZOd\  //daP 


'  NPS  LINAC  /  \dXL 


'Absolute   V     /  deVries 


4 


But 


rd\  (ae& 


^1)  =     d9_  G2 

d"    /Absolute      VdxL     /Mott      d 


therefore 


•i-»(«?L«-^-i 


V  d/i.  y  deVri< 


where  I  ~—  I  ,  T7  .    and  I \  .,  ..  were  calculated  by  programs 


TOPROTON  and  TOPDEUTR  which  are  listed  in  Appendix  E  and  the 
ratio  R  was  the  weighted  mean  given  by  Equation  IV-20. 
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The  statistical  error  in  G,  was  calculated  from 
^/  deVries/\dA  /  Motti 


SR 

c   _  -= ^r_^ 2 

G,  2G, 

d  d 

where  SD   is  calculated  from  Equation  IV-21. 
m 

Several  points  were  taken  in  order  that  the  data  could 

be  checked  for  consistency  on  a  Rosenbluth  plot.   The  results 
are  given  in  Section  V-B. 
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III.   EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

A.  THE  NAVAL  POSTGRADUATE  SCHOOL  LINEAR  ACCELERATOR  (LINAC) 
The  NPS  LINAC  is  a  100  MeV  accelerator  technically  similar 

to  the  Stanford  University  Mark  III  linear  accelerator  which 
is  described  by  Chodorow  et  al_.  [51]  .   A  detailed  description 
of  the  NPS  LINAC  is  given  by  Barnett  and  Cunneen  [52]  and  a 
schematic  floor  plan  is  given  in  Figure  III-l.   The  accelerator 
receives  up  to  22  megawatts  of  peak  power  from  each  of  the 
three  klystrons.   The  machine  is  pulsed  at  a  rate  of  60  pulses/ 
sec,  each  pulse  being  about  1  microsecond  in  duration. 
Electrons  are  injected  at  80  kev  and  can  be  accelerated  to 
energies  between  15  and  100  MeV  with  an  avertige  beam  current 
of  about  5  microamperes.   The  beam  energy  is  selected  by  the 
beam  deflection  magnets  .   The  magnetic  field  of  these  magnets 
is  measured  by  an  NMR  probe.   The  beam  is  focused  down  to  a 
diameter  of  about  2mm  at  the  target  by  the  quadrupole  magnets  . 
The  spread  in  the  beam  energy  is  controlled  by  the  slit 
settings . 

B.  THE  SPECTROMETER  AND  COUNTING  SYSTEM 

The-  spectrometer  is  a  120  ,  16  inch  double-focusing 
magnetic  spectrometer  described  in  detail  by  Oberdier  [44] . 
Ten  detectors  span  a  portion  of  the  focal  plane  of  the  spec- 
trometer, thus  giving  10  counting  channels  for  each  spectro- 
meter setting.   Each  detector  (a  NE102  plastic  scintillator) 
is  7/16  inch  wide  and  1/16  inch  thick.   The  7/16  inch  width 
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allows  n   momentum  spread  Ap/p  of  about  0.3%  for  each  counter. 
Behind  these  ten  detectors  are  two  large  detectors,  both  of 
which  span  the  same  portion  of  the  focal  plane  of  the  spec- 
trometer as  do  the  10  smaller  scintillators.   An  electron 
passing  through  one  of  the  ten  scintillators  produces  a  light 
pulse,  which  is  collected  by  the  photomultiplier  tube,  ampli- 
fied and  transmitted  to  the  control  room  where  it  is  discrimi- 
nated, and  is  then  sent  to  the  scaler  if  both  backing  counters 
simultaneously  registered  outputs.   This  triple  coincidence 
pulse  to  the  scaler  must  also  occur  during  the  time  that  the 
master  gating  circuit  is  pulsing  the  scaler  in  order  that  the 
count  will  register.   The  gating  circuit  prevents  background 
counts  from  registering  between  beam  pulses.   A  block  diagram 
of  the  counting  system  is  given  in  Figure  III-2. 

The  scalers  in  channels  4  through  7  are  capable  of  re- 
cording periodic  pulses  up  to  about  150  M^   while  the  other 

z 
six  channels  of  the  ten  channel  ladder  are  connected  to  scalers 

that  can  record  periodic  pulses  up  to  20  M   .   Only  data  in 

z 
channels  4  through  7  were  analyzed  since  these  could  accumulate 

data  more  rapidly.   The  count  rate  in  these  channels  was  held 
below  60  counts/second  on  all  runs.   Count  rate  correction 
data  were  obtained  and  were  used  to  provide  the  count  rate 
correction  factor  given  in  Equation  (IV-11).  The  correction 
factor  was  linear  out  to  about  100  counts/second,  thus  main- 
taining the  counts  below  60  counts/second  was  more  than 
adequate. 
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The  gas  target  is  covered  in  detail  in  Appendix  D.   Both 
H^  and  d?  were  maintained  at  liquid  nitrogen  temperature  and  at 
about  150  psig  in  order  to  have  a  high  density  of  target 
nuclei. 

C.   METHOD  OF  OBTAINING  DATA 

Consider  first  the  data  which  were  obtained  during  a 
"run"  to  determine  the  ratio  0  /0P  at  a  given  q   and  scat- 
tering angle.   The  same  data  were  taken  for  both  deuterium 
and  hydrogen  and  consisted  of  that  which  was  recorded  auto- 
matically and  a  lesser  amount  that  was  recorded  manually. 

That  which  was  recorded  automatically  consisted  of  any- 
where from  30  to  60  individual  units  of  data,  depending  on 
the  run.   Each  unit  of  data  was  called  a  "point"  and  consisted 
of  the  data  accumulated  at  a  given  spectrometer  setting  while 
a  given  number  of  incident  electrons  passed  through  the  gas 
target.   In  each  run,  the  spectrometer  was  initially  set  such 
that  the  scattered  elastic  spectrum  was  well  above  channel  4. 
The  first  of  the  data  points  was  recorded  there.   Then  the 
energy  of  the  spectrometer  was  lowered  between  0.02  to  0.05 
MeV,  depending  on  the  energy  range  the  spectrometer  was 
sweeping,  and  another  data  point  was  recorded.   This  process 
continued  until  the  spectrometer  setting  was  below  the  peak 
of  channel  7.   The  last  data  point  was  then  recorded. 

After  the  selected  amount  of  beam  charge  had  been  accumu- 
lated at  each  data  point,  the  following  data  were  automatically 
fed  to  a  teletype  machine  that  both  printed  a  hard  copy  and 
punched  a  paper  tape,  a)  the  number  of  counts  recorded  in 
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each  of  the  10  channels,  b)  the  number  of  counts  in  coinci- 
dence in  the  two  backing  counters,  c)  the  time  required  to 
obtain  the  given  charge,  d)  the  incident  energy  of  the 
electrons,  e)  the  capacitance  of  the  capacitor  used  to  in- 
tegrate the  incident  beam,  f)  the  setting  of  the  spectrometer 
in  MeV,  g)  the  maximum  voltage  permitted  on  the  capacitor, 
and  h)  a  number  which,  when  divided  by  30  (the  full  scale 
reading  of  the  integrator) ,  gave  the  fractional  value  of  the 
maximum  voltage  at  which  the  integration  of  the  beam  current 
stopped. 

The  data  which  were  recorded  manually  included  Secondary 
Emission  Monitor  (SEM)  efficiencies  and  the  gas  pressure. 
The  SEM  efficiencies  were  taken  just  prior  to  the  first  data 
point,  at  the  peak  of  each  channel  and  after  the  last  data 
point.   The  gas  pressure  was  recorded  at  the  start  and  finish 
of  each  liquid  nitrogen  "fill"  (as  defined  in  Appendix  C) , 
which  generally  occurred  about  once  at  each  data  point.   The 
leakage  current  in  the  integration  circuits  was  also  recorded 
manually  and  was  used  to  determine  the  ground  loop  errors,  but 
was  not  used  to  calculate  the  cross  section  as  it  was  always 
such  a  small  correction. 

The  above  data  were  always  taken  under  conditions  that 
were  as  identical  as  possible  for  both  the  deuteron  and 
proton.   That  is,  the  experimental  geometries,  the  accele- 
rator beam  current,  the  pressure  of  the  target  gases,  the 
scattered  electron  energy,  the  detection  and  electronics 
systems  and  all  other  possible  variables  were  maintained,  as 
nearly  as  possible,  the  same. 
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2  -2 

Runs  were  made  at  values  of  q   that  went  from  0.02  5  fm 

-2 
to  0.40  fm   .   However,  scattering  from  the  target  windows  or 

_2 

a  contaminant  in  the  gas  masked  the  data  at  0.025  fm   ,  while 

high  background  levels  and  small  cross  sections  made  the  data 

-2 
taken  at  0.40  fm  '  unreliable.   Thus  the  reproducible  data 

-2 

reported  here  covered  a  region  from  0.0  5  to  0.35  fm 

2 
Within  this  range  of  q  ,  forward  angle  scattering  was  prefer- 
entially selected  to  emphasize  the  electric  interaction. 

The  most  significant  change  in  running  conditions  between 
the  deuteron  and  proton  was  the  change  in  beam  energy  made 
so  that  the  scattered  electrons  would  have  the  same  energy  in 
both  cases.   This  insured  that  the  transmission  of  electrons 
through  the  spectrometer  and  the  efficiencies  of  the  detectors 

would  be  the  same  while  taking  D„  and  H„  data.   This  also 

2 

meant  that  the  q   for  hydrogen  was  slightly  different  than 

for  deuterium.   This  difference  was  accounted  for  in  Equation 
11-42. 

The  spectrometer  was  saturated  by  increasing  the  current 
in  the  coils  to  a  setting  corresponding  to  100  MeV  before  data 
were  taken  in  either  H-  or  D~  in  order  to  insure  that  the  same 
hysteresis  loop  could  be  obtained  on  both  spectrums.   The 
energy  defining  slits  were  generally  set  to  give  a  full  width 
at  half  maximum  to  beam  energy  ratio  of  about  0.4%.   The  count 
rate  was  held  below  60  counts  per  second  during  all  runs.   The 
number  of  counts  recorded  varied  from  as  much  as  170,000  to  as 
few  as  10,000  in  a  single  channel.   The  number  of  counts  ob- 
tained depended  primarily  upon  the  size  of  the  cross  section. 
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An.  empty  target  measurement  was  made  at  the  last  spec- 
trometer setting  as  well  as  at  several  others  throughout  the 
energy  range  in  which  data  had  been  taken.   Usually  those  four 

or  five  readings  indicated  a  "flat"  background.   Only  at  the 

2  -2 

lowest  q   of  0.05  F    did  any  structure  show  itself  in  these 

empty  target  background  measurements.   Thus  with  all  but  the 

2  -2 

q  =  0.05  F    data  runs, the  full  target  background  measure- 
ments taken  at  energies  well  above  the  peak  of  each  channel 

were  averaged  and  then  projected  at  those  same  values  through- 

2  -2 

out  each  channel's  elastic  peak.   The  q  =  0.05  F    data  were 

handled  differently  because  there  was  scattering  from  con- 
taminants in  the  gases  and  the  windows  which  complicated  the 
background  spectrum.   In  that  case  empty  target  background  was 
taken  at  each  and  every  spectrometer  setting  of  both  the  hydro- 
gen and  deuterium.   These  empty  target  points  were  then  scaled, 
point  by  point,  to  the  full  target  background  in  order  to  ob- 
tain the  full  target  background  under  each  elastic  peak. 

A  significant  improvement  in  procedure  over  other  work 
which  had  been  done  at  this  facility  is  the  procedure  that 
steps  each  detector  completely  across  each  elastic  peak. 
Formerly,  the  data  from  all  ten  channels  were  used  to  construct 
one  elastic  peak  by  measuring  the  dispersion  between  detectors 
and  the  efficiency  of  all  others  with  respect  to  a  chosen  re- 
ference detector.   This  calculation  of  the  efficiencies  intro- 
duced errors  around  0.5%  into  the  former  method  which  were 
avoided  here. 
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IV..   DATA  ANALYSIS 

A.   GENERAL  METHOD  OF  DATA  REDUCTION 

The  treatment  of  the  experimental  data  (see  Section  III--D) 
is  described  here.   The  first  step  was  to  change  the  form  of 
the  data  from  punched  paper  tape  to  computer  card  data.   This 
was  accomplished  by  a  program  that  required  the  data  in  paper 
tape  form  to  be  read  into  an  on-line  teletype  which  was  con- 
nected to  the  NPS  IBM  360.   This  program  caused  the  data  to 
be  punched  out  on  computer  cards.   Because  of  a  peculiarity 
in  the  data  readout  system  of  the  accelerator,  decimal  points 
were  punched  as  colons.   A  small  program  was  written  which  took 
the  punched  deck  mentioned  above  as  input  data  and  replaced  all 
colons  with  decimal  points  in  a  second  punched  deck.   This  deck 
was  then  interpreted  and  checked  with  the  printed  teletype  out- 
put.  Hydrogen  and  deuterium  data  were  separated  and  were  used 
as  input  to  FORTRAN  programs  T0P#F0RM  where  #  ran  from  channel 
4  through  7.   The  T0P#F0RM  program  is  listed  in  Appendix  D. 
T0P#F0RM  was  written  primarily  to  calculate  the  cross  section 
measured  by  a  particular  channel,  but  it  was  also  used  to 
provide  count  rate  corrected  cata  and  to  normalize  all  points 
to  a  standard  incident  charge,   since  the  automatic  termination 
of  the  integration  at  each  spectrometer  setting  did  not  always 
stop  at  the  same  value.   TOP#FORM  then  prints  out  these  count 
rate  corrected,  charge  normalized  counts  for  each  spectrometer 
setting.   This  part  of  TOP#FORM  was  used  to  plot  these  so- 
called  corrected  counts.   For  that  plotting  run  other  input 

65 


data  which  arc  necessary  to  calculate  a  cross  section,  such 
as  radiation  corrections,  pressure,  etc.  were  given  dummy 
values.   These  count  rate  corrected  and  normalized  counts  for 
all  four  detectors  for  both  the  deuterium  and  the  hydrogen 
were  then  plotted  and  the  position  of  each  peak  and  the  shape 
of  the  scattered  spectrum  was  determined  by  visually  fitting 
a  smooth  curve  through  the  points.   A  typical  plot  of  such 
data  is  given  in  Figures  IV-la  and  IV-lb.   The  difference  in 
energy  between  each  peak  and  the  final  spectrometer  energy  at 
which  data  were  taken  was  then  measured  because  the  radiation 
corrections  depend  upon  Ae,  and  this  value  was  put  into  the 
FORTRAN  program  TOPRADCR  which  is  listed  in  Appendix  E. 
Additionally,  the  peak  energy  (which  is  used  as  the  final 
energy  of  the  electron)  as  well  as  the  pressure  of  the  hydro- 
gen and  deuterium  are  entered  as  input  data  to  TOPRADCR. 
TOPRADCR  has  the  virial  coefficients  of  hydrogen  and  deuterium 
and  the  composition  and  thicknesses   of  the  gas  target  windows 
stored  as  constants.   These  constants  and  the  run  data  men- 
tioned above  permit  TOPRADCR  to  calculate  the  density  of  the 
gas  for  each  channel  and  the  Schwinger,  bremsstrahlung  and 
Landau  radiation  corrections.   The  rationale  behind  the  radi- 
ation correction  calculations  is  given  in  Section  IV-C.   The 
average  background  was  calculated  for  each  channel  by  adding 
all  counts  recorded  prior  to  the  start  of  the  rise  in  each 
channel  and  dividing  by  the  number  of  spectrometer  settings 
which  recorded  these  counts. 
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With  these  calculations  clone,  the  Schv/inger  correction, 
the  product  of  the  Bethe  and  the  Landau  corrections,  the  gas 
density,  the  SEM  efficiency  and  the  background  counts  for 
each  channel  were  typed  on  computer  cards  and  entered  with  the 
appropriate  hydrogen  or  deuterium  deck  as  the  data  for  TOP#FORM 
which  then  calculated  the  elastic  cross  section  and  the  sta- 
tistical error  in  that  cross  section.   The  method  of  calcu- 
lating the  cross  section  is  given  in  Section  IV-B.   The  method 
of  calcultiting  the  statistical  error  is  derived  in  Section  IV-D, 

Four  cross  section  ratios  (one  for  each  channel)  were 

formed  from  the  cross  sections,  and  then  a  weighted  mean  ratio 

was  calculated.   The  electric  form  factor  of  the  neutron, 

Gp  ,  was  then  calculated  by  using  this  weighted  mean  ratio,  R  , 

n 
and  the  appropriate  cross  sections  and  magnetic  corrections 

and  form  factors  as  shown  in  Equations  IV-14  and  IV-15.   The 

values  of  the  proton  electric  form  factors  and  the  proton  and 

deuteron  Mott  cross  sections  and  magnetic  corrections  were 

calculated  from  computer  programs  called  TOPROTON  and  TOPDEUTR 

which  are  listed  in  Appendix  E.   Values  of  D   were  from  data 

in  Appendix  C. 

B.   THE  CALCULATION  OF  ELASTIC  CROSS  SECTIONS 
1.   The  Experimental  Cross  Section 

The  experimental  cross  sections  have  been  calculated 
from  the  expression 

N   K  K     K 

rr  SC    S   b    I     2,-1  ,TT,  t  \ 

ff    =  — t cm  ster  (IV-1) 

Exp    N.   ?7 1  Aft  v 
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where 


N    is  the  number  of  scattered  electrons  counted  in  a 
sc 


particular  channel, 

N.   is  the  number  of  incident  electrons, 
'      in 

2 
77    is  the  number  of  target  atoms  per  cm  , 

Aj\   is  the  solid  angle  defined  by  the  spectrometer 

slits,  and 

JC  •  K      and  ^   are  radiative  corrections  (see  Section 


IV-B) . 


where 


77    is  given  by 

p  t  N 
n     _  a_  (iv-2 

't    A  cos  cp  '  vxv 


2 

p    is  the  density  of  the  target  (gm/cm  ) , 

t    is  the  thickness  of  the  target  (cm) , 

N    is  Avogadro ' s  number, 
a 

A    is  the  atomic  weight  of  the  target  gas. 


and 


<p    is  the  angle  between  the  normal  to  a  plane  solid 
target  and  incident  beam  of  electrons.   In  the  gas  target 
analysis  <p  is  set  equal  to  zero  as  the  angle  cp  is  not  de- 
fined, i.e.,  rotation  of  the  gas  target  still  presents  the 
same  number  of  target  atoms  per  square  centimeter. 

N.   is  determined  by 


in 


where 


N.    =  a  .  ^_  (IV_3) 


Q  =  total  charge  of  incident  electrons 

-9 
e  =  electron  charge  of  1.60210  x  10    coulombs 


70 


c  =  capacitance  of  integrator  in  farads, 

v  =  voltage  of  integrator  in  volts, 

and 

£  =  efficiency  of  SEM. 

Equation  IV-1  is  appropriately  used  to  calculate 

absolute  cross  section,  but  it  was  used  in  this  experiment 

even  though  absolute  cross  sections  were  not  measured.   This 

is  a  correct  procedure  because,  in  this  experiment,  the  ratio 

of  the  cross  sections  of  the  deuteron  and  proton  was  formed 

and  the  systematic  errors  effectively  concelled  as  all  data 

were  taken  under  almost  the  same  conditions  for  both  the 

proton  and  the  deuteron. 

N   ,  the  total  number  of  scattered  electrons,  is 
sc 

obtained  by  integration.   Consider  Figure  (IV-2) ,  where  the 
number  of  detected  electrons  in  a  particular  counting  channel 
is  plotted  against  the  spectrometer  energy. 


and 


Let  E .  be  a  spectrometer  setting 


N.  the  number  of  counts  recorded  in  this  channel 


at  that  spectrometer  setting 
Then 


area  under  scattering  peak 
Nsc  =  (in  Counts  -  MeV) (JV_4) 


where 


$E  =  the  energy  width  of  an  individual  detector. 
This  is  because  each  detector  has  a  finite  dimension 
and  therefore  counts  all  electrons  in  the  spread  Ae.   Thus 
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Number  of 
Electrons 
Counted 


N. 

1 


1+1 


background. 


Spectrometer  Energy  (E) 
Figure  IV-2.  Spectrum  of  Counts  in  a  Single  Detector 


in  obtaining  N   we  have: 
3   sc 


-^tt  dE  where  -5=7  is  the  number  of 

eh  Ctij 


(IV-5) 


counts  per  interval  of  energy.   This  integral  is  numerically 
approximated  by  the  sum 


n-1  h{N.     +  N.  .  ,  )  (E 


N    =  S  , 
SC     1-1 


i+1-  '- 1 
5e 


i+l 


(IV-6) 


The  energy  width  6e  is  obtained  in  terms  of  E.  by  using  the 


dispersion  relation  for  the  spectrometer  optics, 


6e 
E 


fa) 

Dr 


(IV-7) 


where  D  has  been  experimentally  determined  to  be  3.92  [53], 
r   is  the  mean  radius  of  the  spectrometer,  16  inches,  and 
1 W  is  the  width  of  each  detector  (NE102  plastic  scintilla- 
tors) which  is  7/32  inches. 
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Therefore 

6e      7/32" 


E.    3.92(16") 


=  0.0034877  .  (IV-8) 


Putting  these  together,  we  have 

n-1 

N    =  2    MN.-l-N.  J  (E.-E.  ,  )/0. 0034877  E.. 
SC    j^_]_      1   1  +  1     1   l+l  '  1 

The  energy  E.  in  the  denominator  is  the  spectrometer  energy 
setting.   This  setting  is  very  near  the  center  energy  of  the 
band  of  electron  energies  which  the  spectrometer  detects  when 
set  to  the  reading  E..   Hpwever,  the  energy  of  the  electrons 
which  hit  counter  10  will  be  higher  than  any  other  detected 
electrons  and  the  energy  of  the  electrons  which  hit  channel  1 
will  be  the  lowest  in  energy.   In  order  to  take  into  account 
the  fact  that  each  detector  is  recording  electrons  with  dif- 
ferent energies  than  the  spectrometer  setting  E.,  several  runs 
were  made  to  determine  both  the  absolute  calibration  of  the 
spectrometer  and  the  relative  energy  spacing  of  each  of  the 
four  center  channels.   By  comparing  the  calculated  elastic 
peak  energy  with  the  spectrometer  setting  E.,  it  was  deter- 
mined that  the  peak  of  channel  number  5  was  0.17%  higher  than 
the  energy  setting  of  the  spectrometer.   Thus  channel  number 
5  was  chosen  as  the  reference  energy  channel.   By  successively 
lowering  the  spectrometer  setting  from  an  energy  above  the 
elastic  peak  of  channel  number  4  through  channel  number  7,  the 
dispersion  of  the  magnet  was  measured.   This  procedure  was 
done  at  LINAC  energies  of  17  MeV,  30  MeV,  60  MeV,  and  90  MeV. 
The  average  dispersion  of  these  measurements  showed  that  the 
average  energy  of  the  electrons  which  counter  T\    (where  T)    is 
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=  0" 
d-n-     Exp 


4,5,6,  or  7)  detected  was  given  by 

E   =  E.  (1.0  +  0.0039  (77-5))  (IV-9) 

when  the  spectrometer  is  set  at  E . . 
Finally, 

d_a 
Ln. 

(?"  ^[Ni+Ni+l](Ei~Ei+l)/°*0034877Ei  (1-0+0.0039  (ri-5))YsKbKJ. 

/CV    \  (p   t  Na  cm~2\f1    00,    .  _-3   .    ,.    \ 
[— j      [^ -^ ill. 831  x  10    steradians  J 

(IV-10) 
where  the  solid  angle,  Ajx  is  1.831  x  10    steradians. 
2.   Corrections  to  the  Observed  Counts 

a.  Equation  IV-9  would  correctly  measure  an  absolute 
cross  section  if: 

i)  N.,  the  number  of  counts  recorded  on  the 
scalars,  was  the  actual  number  of  electrons  scattered,  i.e., 
if  the  electronic  system  did  not  have  a  finite  resolving 
time.   Even  though  the  scalars  are  capable  of  counting  peri- 
odic pulses  up  to  rates  of  150  MHz,  they  lose  counts  when 
random  (Poisson)  pulses  are  counted  at  high  rates.   At  60 
counts  per  second,  which  was  the  highest  count  rate  which  was 
permitted  during  this  experiment,  about  6%   of  the  counts  were 
lost. 

ii)  N.  consisted  only  of  counts  due  to  the 
scattering  from  the  target  gas  and  was  not  contaminated  by 
background  counts. 

iii)  The  efficiency  of  the  detectors  were  known. 
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b.   The  count  rate  correction  was  determined  experi- 
mentally by  scattering  from  an  aluminum  target  at  count  rates 
up  to  500  counts/sec.   The  number  of  counts  detected  was 
plotted  as  a  function  of  the  count  rate.   Up  to  about  100 
counts/second  the  number  of  counts  recorded  decreased  linearly 
with  rate.   A  linear  fit  was  made  to  the  data  in  this  region, 

and  gave   ._ .  .  ,         .  ,»    N  (recorded)      ,_„  ,  ..  , 

N (count  rate  corrected)  =  iIoTo0097N/t     (1V-11) 

where  N/t  is  the  average  count  rate  determined  by  dividing  the 
recorded  number  of  counts  by  the  time  of  the  integration. 

c.  The  background  was  determined  by  counts  obtained 
at  energies  greater  than  the  elastic  peak  of  either  H2  or  D~, 
whichever  was  the  target  at  that  time,  for  each  of  the  four 
detectors 0   The  average,  number  of  counts  for  each  channel,  in 
this  region,  was  then  determined  and  was  used  as  the  background, 
During  each  run  an  empty  target  background  was  also  taken  to 
insure  that  there  was  no  rise  or  fall  in  this  extrapolated 
average  background  under  the  elastic  peak  caused  by  some  source 
other  than  the  H-  or  D^. 

d.  The  efficiency  of  each  detector  was  not  deter- 
mined because  the  object  of  this  experiment  was  to  obtain  the 
ratio  of  the  deuteron  cross  section  to  the  hydrogen  cross 
section.   Because  this  ratio  was  determined  for  each  of  the 
four  detectors  separately,  the  efficiency  of  each  detector  can- 
celled and  therefore  is  not  included  in  the  determination  of 
the  cross  section. 

e.  An  electrometer  was  used  to  measure  the  inte- 
grated incident  beam  current  for  each  data  point.  The  cir- 
cuitry associated  with  the  electrometer  was  set  to  stop  the 
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integration  and  the  counting  system  whenever  a  predetermined 

voltage  across  a  capacitor  charged  by  the  current  from  a  SEM 

was  reached.   The  electrometer  did  not  always  stop  at  exactly 

the  same  voltage  but  varied  sometimes  by  as  much  as  a  percent 

or  two.   Therefore  both  the  hydrogen  and  the  deuterium  counts 

were  normalized  to  a  standard  voltage  (V„  £       )  for  each 

J     Reference 

run.   Thus,  if  V.  ,   _  was  the  actual  voltage  at  the  end  of 

Actual 

the  integration,  then  the  number  of  counts  recorded  was  cor- 
rected by 

N  (corrected  for  voltage)  =  Na  .   ,   Reference  .  (iv-12) 

Actual  VActual 

The  cross  section,  after  correcting  for  count  rate, 
background  and  voltage  cutoff  is 

V 


.n-1 


BKG 


^(l-0,00097Ratei) 


ref 

Ni+1   Vact BKG 

(l-0.00097Ratei=j) 


d"      /CV 


\e£    4 
SEM 


ptNacm   \  (.    __.    -_-3       , .    ^ 
z. j    ^i#831  x  10    steradiansy 


(E.-E.  .)K  K.K 

1 1+1 sjb_l_  ,  fTV-13 

0.0034877  (l.O-i-0039  (??-5)  )E.  ; 


C.   RADIATION  CORRECTION  CALCULATIONS 

The  proton  cross  section  was  derived  by  Rosenbluth  [26] 
with  the  assumption  that  a  single  virtual  photon  is  exchanged 
in  the  scattering  process  as  shown  in  Figure  IV-3. 

The  deuteron  cross  section,  as  derived  by  Jankus  [13]  in 
the  impulse  approximation,  also  assumes  that  one  virtual  pho- 
ton is  exchanged  as  shown  in  Figure  IV-4. 
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electron 


Feynman  Diagram  of  Rosenbluth  Calculation 
Figure  IV-3 


Feynman  Diagram  of  Impulse  Approximation 
Figure  IV-4 

Unfortunately,  when  one  does  an  electron  scattering 
experiment  these  processes  are  not  the  only  ones  which  occur, 
and  the  experimenter  must  take  into  account  competing  pro- 
cesses in  order  to  be  able  to  compare  experiment  and  theory. 
In  fact  the  experimental  scattered  energy  spectrum  which  is 
measured  contains  the  effects  of  three  such  competing  pro- 
cesses which  must  be  accounted  for  before  one  can  obtain  the 
cross  sections  for  processes  as  shown  in  Figures  IV-3  or  IV-4 
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These  competing  processes  are: 

1)  Photon  emission  and  absorption  while  the  electron  l's 
in  the  field  of  the  scattering  nucleus  of  interest..  The  .cor- 
rection to  the  cross  section  due  to  this  process  was  first 
calculated  by  Schwinger  [54]  and  is  commonly  called  the 
Schwinger  Correction. 

2)  Bremsstrahlung  emission  in  the  target  before  or  after 
the  nuclear  scattering  of  interest. 

3)  Excitation  and  ionization  of  material  in  the  target. 
Each  of  the  above  processes  causes  the  energy  of  a 

scattered  delectron  to  be  less  than  it  would  be  if  the  only 
process  which  occurred  was  the  one  photon  exchange.   Thus,  a 
radiative  tail  appears  on  the  low  energy  side  of  the  elastic 
peak,  as  shown  in  Figure  IV- 5. 


Energy  of  Scattered  Electrons  "• 

Elastic  Peak  with  Radiation  Tail 
Figure  IV-5 
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To  obtain  a  cross  section  for  the  nuclear  scattering 
event  of  interest,  the  experimenter  would  have  to  continue 
down  the  radiative  tail  to  zero  energy  in  order  to  capture 
those  scattered  electrons  which  have  some  probability  of 
losing  all  their  energy  in  the  three  processes  listed  above. 
In  practice,  the  measurement  is  stopped  at  Ae  below  the 
elastic  peak  (see  Figure  IV-5) .   The  value  of  Ae  is  generally 
of  the  order  of  at  least  four  half-widths  (T/2) [55] .   Thus 
fewer  electrons  are  detected  and  the  experimental  cross 
section  is  less  than  it  should  be.   Since  each  of  the  three 
processes  mentioned  above  contributes  to  the  radiative  tail, 
corrections  must  be  made  to  account  for  the  undetected  elec- 
trons with  energy  less  than  Ef  -  Ae. 

1 .   Schwinger  Correction 

The  largest  of  the  corrections  is  that  which  accounts 
for  photon  emission  and  absorption  while  the  electron  is  in 
the  field  of  that  nucleus  which  causes  the  large  angle  scat- 
tering which  is  of  interest.   This  correction  to  the  cross 
section  was  generally  about  15%  in  this  experiment. 

The  Schwinger  correction  accounts  for  many  inter- 
actions which  may  occur  in  addition  to  the  one  virtual  photon 
exchange  shown  in  Figures  IV- 3  and  IV-4.   Several  of  these 
possible  interactions  are  shown  in  Figure  IV-6. 

In  processes  (a)  and  (b)  real  photons  are  emitted 
before  and  after  the  nuclear  scattering.   Processes  (c) ,  (d) , 
and  (e)  are  called  vertex  renormalization  terms  and  (f)  is  the 
vacuum  polarization  term.   Improved  radiative  corrections  to 
electron-proton  scattering  have  been  calculated  by  Tsai  [56] 
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Feynman  Diagrams  of  Schwinger  Correction  Interactions 
Figure  IV-6 

and  by  Meister  and  Yennie  [57]  and  include  the  effects  of 
proton  recoil.   Numerical  evaluation  of  the  Schwinger  cor- 
rections in  this  work  used  the  results  given  in  the  review 
article  by  Maximon  [58] . 
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where  E.  and  Ef  are  the  electron  energies  prior  to  and  after 

scattering,  ol   is  the  fine  structure  constant,  m  is  the  elec- 

2 
tron  mass,  q   is  the  four-momentum  transfer,  and 
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l2(x)  =  -  $X  i<*Mi-t).  at 

O 

is  the  Spence  function. 

2 .   Bremsstrahlung  Correction 

When  passing  through  the  target  material,  the  elec- 
trons also  radiate  away  energy  through  radiative  collisions 
with  atomic  electrons  and  nuclei  other  than  the  target  nuc- 
leus of  interest.   The  energy  loss  due  to  bremsstrahlung  has 
been  studied  by  Bethe  and  Heitler  [59]/  who  calculated  the 
cross  section  $ (E  )  and  the  average  energy  loss  due  to  this 
process.   By  approximating  this  cross  section  Heitler  [60]  and 
Bethe  and  Ashkin  [61]  obtain: 

dE 
_  _o  _  „  „ 

dx   "   o  o 

where  X   is  the  radiation  length,  which  is  the  distance  over 
which  the  electron  has  its  energy  reduced  by  a  factor  1/e. 
In  this  experiment,  the  expression  given  by  Butcher  and 
Messel  [62]  has  been  used  to  calculate  the  radiation  length 

V 


X 


°    40WQ  r^z(z+0.8)p  0n(183z  1//3) 

where  N   is  Avogadro's  number  and  r   is  the  classical  elec- 
o  e 

3 
tron  radius,  p  is  the  density  in  gm/cm   and  A  is  in  grams/mole. 

The  target,  which  consisted  of  the  hydrogen  or  deu- 
terium gas,  was  in  a  pressure  vessel  which  had  0.001  inch  steel 
windows  for  the  incident  and  scattered  beam.   For  the  brems- 
strahlung correction,  and  also  the  Landau  corrections,  the 
target  was  treated  as  a  homogeneous  molecular  compound  of 
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steel  and  either  1I„  or  D„  which  consisted  roughly  of  Fe, 

H~. ,  or  Fe.  D,.r,  where  the  "molecular  formula"  was  calcu- 
26        1   26 

lated  for  each  gas  pressure.   This  formula  accounts  for  the 
1  mil  windows  and  the  2  inches  of  target  gas  at  150  psig  and 
77  K.   Following  Butcher  and  Messel  the  radiation  length  for 
such  a  compound  target  was  calculated  from 


1 
X 


In  the  actual  FORTRAN  programs  the  X  were  multiplied 

2 
by  the  average  density  of  this  compound  to  give  X   in  gm/cm  . 

Bethe  and  Heitler  [59]  also  give  the  probability  that 

an  electron  traversing  t  radiation  lengths  will  have  an  energy 

loss  less  than  a  times  its  initial  energy  (aside  from  the 

ionization  loss  and  the  scattering  from  the  nucleus  of 

interest) .   That  probability  is 


a(t/Zn  2) 


c  (a)   r(i+  t/&i  2)  * 

Hofstadter  has  integrated  this  and  calculated  the  cor- 
rection factor  to  the  observed  cross  section;  this  brems- 
strahlung  correction  factor  may  be  written  as  [6] 

\  =  exp  (x~!r2  )log  {7k]   ' 

o 

where  T  is  the  thickness  of  the  target,  including  the  two 

2 
steel  0.001  inch  windows  as  well  as  the  gas,  in  gm/cm  ,  and 

E  is  the  geometric  mean  electron  energy  while  going  through 

the  target.   PL,  was  about  1.02  in  this  experiment. 
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3.   Ionization  Correction 

As  the  electrons  in  the  beam  pass  through  the  target 
they  also  lose  energy  by  excitation  and  ionization  in  the  gas 
and  in  the  gas  target  windows.   Landau  [63]  first  calculated 
both  the  average  energy  a  monoenergetic  beam  would  lose,  and 
the  spread,  or  straggling,  of  the  energy  distribution  of  the 
beam  after  passing  through  a  thin  target.   From  this  distri- 
bution it  is  possible  to  calculate  the  fraction  of  electrons 
which  are  not  detected  because  data  is  taken  only  down  to 
E  -Ae.   Figure  IV-7  shows  a  theoretical  spectrum  of  monoener- 
getic incident  electrons  after  passing  through  a  layer  of 

material.   T  is  the  full-width  at  half -maximum.   0   v  is  the 

prob 

most  probable  energy  loss  due  to  ionization  and  E     is  the 
incident  energy,  and  N(E)  is  the  number  of  electrons  detected 
at  energy  E. 


N(E) 


Eo  Ei  El 


INC 


Effect  of  Ionization  on  Spectrum  of 
Monoenergetic  Incident  Beam 

Figure  IV-7 
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Following  Landau's  [63]  method,  but  using  the  general 

notation  of  Breuer  [64],  one  can  calculate  the  fraction,  Y(Ae) 

of  the  beam  one  does  not  detect  (because  of  ionization  losses) 

if  data  taking  stops  Ae  below  the  peak  in  Figure  IV-7  as 

E„  E. 

r-   2  n   mc 

Y(AE)  =  ^   N(E)dE  /  ^     N(E)dE  . 

o  o 

Landau  [63]  gave  a  graph  of  ^  (Ae)  but  he  also  derived 

an  approximation  for  Y(Ae),  as  have  others  [65],  which  has 

been  used  in  this  work: 

where  the  value  of  Ui  as  provided  by  B6rsch-Supan  [66]  is 
co  =  X(l-  -r+T~  ) 


where  A  =  y    -    1  =  -0.4228..., 
y   is  Eulers  constant, 

Ae 

X  =  -*-  -  0.05  (is  an  approximation  to  the  X  as  defined 
by  Landau  and  Bergstrom  [67], 

£  rj.z. 

?  =  0.154  v  J-^x   T    MeV,  and 

'l  1 

2 
T  =  target  thickness  is  gm/cm  . 

77.  is  the  number  of  atoms  of  element  (i)  per  molecule 

of  the  target  which  here  was  considered  to  be  a  homogenous 

mixture  of  the  0.001  steel  windows  and  the  gas  in  the  target 

vessel,  and  z .    and  A.  are  the  atomic  number  and  atomic  weight 

of  element  (i) . 

The  ionization  correction  is  given  by 

K.  =  (1-  Y(AE))"1 
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In  this  experiment  K.  V/as  less  than  1.01,  i.e.,  the 
Landau  ionization  correction  was  less  than  one  percent. 

D.   ERROR  ANALYSIS 

1 .   Systematic  Errors 

Because  this  experiment  measured  the  ratio  of  the 
deuteron  to  the  proton  cross  section,  many  systematic  errors 
in  each  of  these  cross  sections,  such  as  solid  angle,  scat- 
tering angle,  SEM  efficiency,  target  thickness,  and  counting 
rate  corrections  were  effectively  cancelled  when  the  ratio  of 
the  two  cross  sections  were  formed.   This  follows  because  the 
data  for  both  cross  sections  were  taken  during  one  continuous 
LINAC  running  period  v/ith  the  geometries  the  same  and  the  beam 
characteristics,  excluding  the  incident  energy,  as  nearly 
identical  as  possible. 

A  systematic  error  which  does  not  cancel  in  this 
ratio  experiment  is  that  which  can  arise  because  of  possible 
impurities  in  the  hydrogen  and  deuterium  gases.   The  supplier 
of  the  target  gases  advertized  99.5%  pure  deuterium  and 
99.995%  pure  hydrogen.   However,  after  closer  investigation 
it  was  discovered  that  their  analytic  equipment  can  only  de- 
tect 0.3%  or  more  contamination  of  D?  in  H2  or  H2  in  D„.   All 
the  data  taken  in  this  experiment  used  only  one  bottle  of 
deuterium  and  one  of  hydrogen,  thus  all  data  were  taken  with 
the  same  amount  of  possible  contaminant.   If  there  were  a 
difference  in  the  contamination  of  H2  or  D2  this  difference 


would  have  caused  a  non  zero  intercept  of  the  plot  of  Gp 

2  n 

versus  q  .   For  example,  a  0 . 3%  contamination  of  D2  and  a 
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0.0%  contamination  of  H_  would  have  lowered  the  intercept 
by  about  0.0015. 

Whenever  empty  target  background  was  compared  to  the 
background  above  the  peak,  the  deuterium  background  always 

had  the  same  spectrum  distribution  in  both  the  empty  and 

2  2 

full  cases.   However,  at  lower  q   (particularly  q   =  0.05) 

the  hydrogen  spectrum  above  the  hydrogen  peak  clearly  showed 

a  contaminant  in  the  gas  which,  if  it  were  a  low  atomic 

number  contaminant  would  have  been  a  few  percent,  but  if  it 

were  a  high  atomic  number  contaminant  it  would  have  been  less 

than  a  few  tenths  of  a  percent.   The  closeness  of  the  inter- 

2 

cept  of  G    versus  q   to  zero  (we  report  0.0000  +  0.0012  for 

n 
LF  model  15)  indicates  that  the  H?  and  D2  contaminants  were, 

in  fact,  almost  equal.   The  following  figure  shows  the  dif- 
ference between  the  H2  and  D2  spectra.   The  elastic  peaks  in 
these  cases  would  be  at  about  1000  counts. 
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2 .   Random  Errors 

The  random  errors  in  this  ratio  experiment,  besides 
the  statistical  (Poisson)  counting  errors,  came  from  errors 
in  determining  the  density  of  the  gas  in  the  target,  vari- 
ations in  the  SEM  efficiency  during  a  run,  energy  shifts  of 
the  incident  beam,  leakage  currents  in  the  beam  current 
integrators,  and  very  small  errors  in  the  radiation  cor- 
rections resulting  from  errors  in  determining  Ae. 

The  statistical  error  for  each  cross  section,  due  to 
the  Poisson  distribution  of  the  number  of  electrons  scat- 
tered for  a  given  number  of  incident  electrons,  was  calcu- 
lated by  a  subroutine  in  the  same  program  that  calculated 
the  cross  section  and  included  the  statistical  errors  made 
in  determining  the  background  subtraction.   To  determine  the 
fractional  error  in  the  cross  section  caused  by  these  sta- 
tistical fluctuations,  the  cross  section,  X,  may  be  written 

as  M 

X  =  Constant  •  2   (N,  -B) 
k=l   K 

where 

N,  =  Total  counts  (count  rate  and  voltage  corrected) 
obtained  at  the  k    spectrometer  setting, 
B  =  Average  background  to  be  subtracted  (determined  at 
energies  above  the  peak)  and  assumed  to  be  the  same 
for  each  spectrometer  setting,  and 
M  =  The  number  of  settings  of  the  spectrometer. 

M 

Let  N  =  2   n,  be  the  total  number  of  counts,  then: 
k=l  K 

X  =  Constant  •  (N-M'B) . 
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The  variance  of  the  cross  section  is 

^   2 

C©/   "B 


CT2 .  (&Y  a2  +  (^y  a 

x  \bN/         N 

which   gives 


tf2  =    (Constant)  2  (N+M2^2) 

2  2 

where  ct   is  the  variance  in  the  average  background  and  cr   =  N. 

Let  B   =  ?_.  B.  (where  B.  =  background  at  i  ^  setting)  and 
J  is  the  number  of  spectrometer  settings  above  the  peak  at 

which  background  was  taken. 

2b  J  J 

Since  B  =  -ji  then  og  =  -£  S  eg  =  -t  I  6,=  i  BT  -  §. 

^^  2  2  2 B 

cr  =  (Constant)  [N+M   ~  ] 

and  the  fractional  standard  deviation  of  the  cross  section  is 


cr     r        2     .    >h 

x  _  JN+M  B/J]  2 

X     [N-MB] 
which  can  also  be  expressed  as 

2 

x 1_    MB  (1+  M/J) 

X2  =  No       N2 
o 

where  N   =  N-MB  and  the  first  term  gives  the  variance  due  to 
o  ^ 

the  true  counts  and  the  second  gives  the  variance  resulting 
from  the  background. 

The  errors  in  the  measured  cross  sections  caused  by- 
errors  in  determining  the  density  of  the  gas  in  the  target 
varied  slightly  from  run  to  run,  but  were  always  less  than 
about  0.1%.   The  pressure  dropped  while  liquid  nitrogen  was 
being  pumped  and  then  rose  again  and  stayed  steady  until  the 
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next  nitrogen  fill,  but  the  variations  were  almost  identical 
in  the  hydrogen  and  the  deuterium,  and  were  about  0.3%  of  the 
measured  pressure.   In  calculations,  the  average  pressure  at 
the  start  of  the  nitrogen  fill  or  at  the  completion  of  the  fill 
was  the  pressure  actually  used,  depending  on  which  was  more 
constant.   In  calculating  the  gas  density,  of  course,  if  the 
pressure  at  the  start  (completion)  of  the  fill  was  used  for 
hydrogen,  then  the  pressure  at  the  start  (completion)  of  the 
fill  was  also  used  for  deuterium. 

Generally,  six  SEM  efficiencies  were  measured  during 
a  data  run  on  each  gas;  one  was  taken  at  the  start  of  the  run, 
one  at  each  of  the  four  peaks,  and  one  at  the  completion  of 
data  taking.   At  each  point,  four  or  more  measurements  were 
made,  and  the  fractional  standard  deviation  of  a  single  mea- 
surement was  usually  less  than  0.2%,  giving  a  fractional 
standard  deviation  of  the  mean  of  the  SEM  readings  of  less 
than  0.08%  for  the  six  readings. 

Because  of  problems  such  as  variations  in  line  voltage, 
the  beam  energy  was  not  constant  and  thus  the  scattered  energy 
also  would  change  slightly  from  time  to  time  throughout  a  run. 
The  result  of  these  changes  was  a  shift  of  counts  from  one 
channel  to  another  of  the  counting  ladder.   This  was  most 
noticeable  when  data  were  being  taken  along  the  steep  slopes 
of  elastic  peaks;  often  a  low  count  in  one  channel  would  be 
accompanied  by  a  high  count  in  an  adjacent  channel.   The 
obvious  method  of  reducing  these  shifts  was  to  close  the 
energy  defining  slits  more,  but  this  increased  the  background 
and  thus  the  ratio  of  real  counts  to  background  became 
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less  -  giving  larger  errors  on  the  measured  cross  sections. 
Nevertheless,  the  experimenter  could  not  always  unambiguously 
estimate  how  many  counts  were  added  or  subtracted  from  a  given 
elastic  peak  due  to  energy  shifts. 

An  examination  of  many  runs  shows  that  there  was  no 
apparent  energy  shift  problem  in  those  runs  where  the  cross 
section  was  large  and  the  slits  could  be  closed  and  still 
maintain  a  good  peak  to  background  ratio.   In  other  runs 
where  the  slits  were  opened  wider,  it  appeared  that  a  varia- 
tion of  as  much  as  a  percent  could  have  occurred  at  several 
points  in  an  elastic  peak.   Adjacent  data  points  and  general 
elastic  peak  shapes  were  used  to  smooth  out  these  large  fluc- 
tuations.  It  is  estimated  that  in  the  poorer  runs  these 
energy  shifts  could  have  caused  an  error  in  the  average  mea- 
sured cross  section  ratio  of  as  much  as  0.2%. 

There  were  slight  leakage  currents  in  the  circuits 
used  to  measure  and  integrate  the  beam  current.   These  cur- 
rents were  reduced  by  improving  the  insulation  of  various 
components  and  by  cleaning  all  sensitive  elements.   The 
leakage  currents  on  the  SEM  and  the  faraday  cup,  along  with 
all  the  circuits  used  to  accumulate  the  currents,  were  mea- 
sured and  found  to  be  almost  trivial.   These  measurements 
were  repeated  during  each  run  and  remained  very  constant.   If 
the  electron  beam  was  perfectly  constant  during  all  parts  of 
a  given  run,  the  leakage  current  would  introduce  no  errors. 
Usually  the  electron  current  was  held  constant  to  within  20%; 
this  would  cause  an  error  in  the  ratio  of  the  cross  sections 
of  about  0.02%. 
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Finally,  there  is  the  error  in  measuring  the  exact 
position  of  the  peak  of  the  elastic  spectrum  from  the  plotted 
data.   The  amount  of  this  error  varies  with  such  things  as 
beam  stability,  energy,  and  the  widths  of  the  slits  but 
generally  was  about  0.01  MeV.   This  error  causes  an  error 
in  the  measurement  of  Ae  which  in  turn  causes  an  error  in  the 
Schwinger,  bremsstrahlung,  and  ionization  corrections.   The 
error  in  the  measured  cross  section  due  to  this  Ae  error  for 
all  three  radiation  corrections  was  of  the  order  of  0.05%. 
The  error  in  Ae  also  causes  an  uncertainty  in  the  final  energy 
value  at  which  one  terminates  the  integration  of  the  elastic 
spectrum.   This  error  caused  about  a  0.03%  error  in  the  mea- 
sured cross  section. 

Combining  all  the  above  errors  except  the  stat.i  stical 
errors  and  the  energy  shift  errors  (since  these  two  vary  so 
much  from  run  to  run) ,  it  is  seen  that  the  total  fractional 
error  in  the  measured  cross  section  is  about  0.14%.   In 
addition  to  this,  the  statistical  errors  in  the  cross  sections 
varied  from  an  average  of  less  than  about  0.2%  to  a  maximum  of 
about  0.4%  for  a  few  runs,  while  the  energy  shift  errors  varied 
from  apparently  zero  through  an  average  of  less  than  0.1%  with 
about  4  runs  possibly  having  errors  as  large  as  0.2%.   A  sum- 
mary of  the  errors  discussed  above  is  given  in  Table  IV-1. 

The  average  error  in  the  cross  section  ratio  due  to 
all  random  errors  except  the  statistical  error,  S,    _.,    , 
has  been  calculated  from  the  values  of  the  average  cross 
section  errors  in  Table  IV-1  to  be  0.24%.   The  total  error  of 
the  weighted  mean  ratio  R  was  then  formed  from 
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2         2         2 
STotal  =  SStat.  +  SAvg. Others  *  (IV-1) 


TABLE  IV-1 

Random  Errors  in  Cross  Section  Data 

Errors  Average  Values  of  Errors 
Statistical  0.3% 

Gas  Density  0.1% 

SEM  Efficiency  0.08% 

Energy  Shifts  0.10% 

Leakage  Currents  0.014% 

Ae  Measurement  (Radiation  Corrections)  0.05% 

Ae  Measurement  (Integration)  0.03% 

Two  significant  sources  of  error  were  avoided  in  this 
work  that  were  present  in  previous  work  [2] .   The  use  of  a  gas 
target  in  place  of  a  polyethylene  target  and  the  independent 
determination  of  the  cross  section  ratios  by  each  channel  per- 
mitted these  improvements.   The  gas  target  scatters  electrons 
from  the  target  nuclei  all  along  the  beam  path,  but  only  those 
nuclei  which  satisfy  the  scattering  geometry,  i.e.,  those 
which  will  scatter  and  recoil  through  both  a  collimator  which 
is  about  6  inches  from  the  gas  target  and  the  spectrometer 
slits  which  are  16  inches  from  the  target,  can  enter  the  spec- 
trometer.  Thus  if  the  experimenter  first  scatters  from  hydro- 
gen, maintains  the  exact  geometry,  and  then  fills  the  target 
with  deuterium  and  scatters  from  it,  he  does  not  have  to 
calculate  the  target  thickness  in  order  to  obtain  a  perfectly 
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accurate  ratio  of  cross  sections  -  if  target  thickness  were 
the  only  possible  error  in  the  cross  section  determination. 
This  is  not  the  case  if  he  used  a  solid  target  which  is 
thinner  than  the  collimated  entrance  width  to  the  spectro- 
meter.  Stewart  [2],    who  used  a  thin  polyethylene  target  to 
measure  the  deuteron  and  proton  cross  sections,  estimated 
the  target  thickness  error  as  0.5%.   A  second  significant 
error  which  was  avoided  in  this  work  was  the  error  experi- 
enced in  determining  the  relative  efficiencies  of  each  of 
the  ten  channels.   This  error  was  quite  informally  estab- 
lished to  be  0.5%  in  Ref.  2.   In  this  present  work  each  de- 
tector scanned  the  complete  elastic  spectra  and  the  ratio 
of  the  deuteron  cross  section  to  the  proton  cross  section 
was  formed  separately  for  the  individual  detectors.   This 
eliminated  any  detector  efficiency  errors  which  would  have 
arisen  had  the  method  of  Ref.  2  been  used. 

Having  established  what  errors  were  made  in  the 
measured  cross  sections,  one  can  determine  how  these  errors 

propagated  into  the  value  of  G   .   Let  the  standard  deviation 

n 
of  a  cross  section  be  denoted  by  S   and  S  ,  for  the  proton  and 

deuteron.   Let  R  =  O   /ap  where 

cr   =  experimental  cross  section  of  the  deuteron  and 

0"^  =  experimental  cross  section  of  the  proton. 

G_,   is  determined  from 

,\  \ 

ge  -  Uri-1-0  H  (IV_14) 

n   V  Ep   C       I       p 
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and  G   /G    is  determined  from 
Ed 

aP     Kd    GE  (q^ 
R_T4ott   ]W_  _^_|_   .  (IV_15) 

«Mott   *£     GE  «£> 
P 

Solving  first  for  the  fractional  error  in  R  due  to 
the  errors  in  both  cross  sections  we  have  [42] : 


R     (f  )   .  (°P) 


where  S   is  the  standard  deviation  of  the  ratio  R. 

R  GF 

Ed 
The  error  in  the  ratio  of  -■  ■  ■   is  determined  from 

GE 

Equation  1V-15,  thus  the  standard  fractional  deviation  of  the 

ratio  G_  /G_  ,  denoted  by  S"    ,         /is 
Ed   Ep  S/  Ep 

S*  ,    -%£.  (IV-17) 

d   P 

From  Equation    IV-14  we   determine    the   variance    in   G„ 

n 
to  be:  v  2  /G       (a2\ 

and  using  Equation  IV-15,  the  standard  deviation  in  the 
neutron  electric  form  factor  is 

aP  Kd  \  k  °E  1%) 

_  L  o  I  1  Mott  _JB£  P  1         /TV  1q 

SG   -  *  SR   R  ad       p   J  D  *        (IV~19 

En        \   Wt  Km  /  C 


Up  to  this  point  in  the  error  analysis  only  the 

errors  in  one  channel  have  been  discussed,  i.e.,  S   is  the 

R 
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error  in  R  due  to  errors  in  the  dcuteron  and  proton  cross 

sections  for  any  given  channel.   Letting  R.  be  the  ratio 

measured  by  channel  i  and  S  .  its  standard  deviation,  then 

Ri 

the  weighted  mean  of  all  four  channels,  R  ,  for  any  given 

run  is  obtained  by  calculating 

R. 
1 
7 
2 

Rm=   f 4  /_1^ 
and  the  standard  deviation  of  this  weighted  mean  ratio  is 


m    -(^ 

R   and  S    are  the  experimental  values  of  the  ratio 
m 
and  of  the  standard  deviation  of  that  ratio  quoted  in 

Table  V-I  for  each  run  and  are  the  values  used  to  determine 

all  derived  quantities  for  each  run  in  this  work,  e.g., 

2 

G„    ,    S_      ,    G,    and   all    others. 

n  E 

n 


E.   REJECTION  OF  DATA 

Twenty  one  runs  are  reported  in  this  work.   In  each  run 
the  cross  section  ratio,  R,  was  measured  for  each  of  the  four 
channels,  giving  a  total  of  84  measurements  of  R.   However, 
nine  of  these  84  measurements  have  not  been  included  in  the 
data  analysis  because  of  clear  indicators  that  each  had  large 
systematic  errors. 
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The  data  given  in  Table  V-I  for  q  =  0.05  at  0  =  60° 
does  not  include  the  data  of  channel  5  since  the  "contamina- 
tion" background  subtaction  method  discussed  in  Section  JII-B 

2 
for  the  q   -  0.05  data  did  not  work  well  for  that  channel. 

Energy  shifts  caused  clear  irregularities  in  four  runs. 

For  this  reason  the  following  were  not  included  in  the  data 

2 
reported  in  Table  V-I:  channel  5  in  the  second  q   =0.10  run, 

2 

channel  G  in  the  first  q   =0.2  run,  channel  7  in  the  second 

q   =0.2,  0  =  90   run,  and  channels  4,  5,  and  6  in  the  second 
q2  =  0.2,  0  =  60°  run. 

The  inability  to  establish  reliable  backgrounds  caused 
the  rejection  of  the  data  of  channel  4  in  both  q   =0.35 
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V .   EXPERIMENTAL  RESULTS 

A.   NEUTRON  CHARGE  FORM  FACTORS 

The  data  from  21  runs  are  reported  here.   Each  run  con- 
sisted of  four  measurements  of  R,  the  ratio  of  the  elastic 
electron-deuteron  cross  section  to  the  elastic  electron- 
proton  cross  section.   Of  these  84  measurements  of  R,  9  are 
not  included  in  the  reported  data  for  reasons  given  in 

Section  IV-E. 

2 
Table  V-I  gives  the  q   transferred  to  the  deuteron,  the 

scattering  angle  9,  the  incident  energies  of  the  electron 

incident  on  the  proton,  E.  ,  and  deuteron,  E.  ,  the  experi- 

1p  xd 

mental  ratio  R   and  the  standard  deviation  of  the  ratio  R  as 

m 

determined  by  the  total  error  (defined  in  Section  IV-D)  and 
the  statistical  (Poisson)  error  for  that  run.   Also  given  are 
quantities  which  were  determined  by  calculations.   These 
quantities  are  explained  in  Section  II.   Using  the  experi- 
mental results  and  the  calculated  quantities  given  in  Table 
V-I,  the  ratio  of  the  deuteron  electric  form  factor  to  the 
proton  electric  form  factor  has  been  derived  and  is  listed 
in  the  right  hand  column. 

The  neutron  electric  form  factor  G   was  obtained  by 

hn 
using  Equation  11-43  in  the  following  form 

crP(G)    Kd\  M   GE     {4]    , 

mr  \      )         p      f*     (     _1_ 


n 


°X  j  )\^\  Ce 


97 


-    1 


Vq* 


The  deuteron  charge  form  factor  C   listed  in  Appendix  C  and 

til 
2 
the  G   /G    and  G   (q,)  values  given  in  Table  V-I  were  used 

d    P        P  • 


to  solve  for  G   .   The  value  of  C   depends  upon  the  deuteron 

n  2 

wave  function  as  well  as  q  .   Three  sets  of  the  21  values  of 


G   have  been  calculated  for  Lomon-Feshbach  models  #1,  #5, 
n 

and  #15. 

The  best  estimate  of  the  error  in  each  of  these  21  values 

of  G   ,  for  each  LF  model,  is  the  error  obtained  using  Sm    . 
xj  j.o  tax 

from  Equation  IV-1.   The  best  estimate  of  the  value  of 

G_   reported  for  each  of  the  three  models  was  obtained 

-,    /.      jti 
dq     n 

from  these  sets  of  G„   weighted  by  S„  ,  , .   For  comparison, 
rij  i  o  t  ax 

7T   G   was  also  calculated  from  the  G„   data  weighted  by  its 

n   £■  E  111 

dq     n  n 

statistical  error  only,  and  also  with  equal  weighting  given 

to  each  GE   to  provide  an  insight  into  the  sensitivity  of  the 

n 
result  to  weighting  of  the  data. 

After  calculating  each  G„   and  its  Sm  .  ,  an  average  G„ 
3        E  Total  J    E 

P     n  n 

and  its  Srp  ,  ,  at  each  q  was  calculated.   A  plot  of  the 

2 

average  Gp   (with  its  S™  ±.a->)    versus  q'  is  provided  for  each 

n 

of  the  three  models.   For  LF  model  15  the  value  of  each  G„ 

E 
n 

and  its  S_  .  ,  is  also  given. 

For  the  21  rums  reported,  all  values  of  Gp   were  calcu- 

n 
lated  from  the  FORTRAN  programs  listed  in  Appendix  E.   In 

addition,  in  order  to  check  these  programs  in  15  runs,  the 

values  of  G„   were  also  calculated  on  a  Wang  Model  700  cal- 

En 
culator  using  programs  not  written  by  the  author.   These 

Wang  calculations  used  the  same  raw  data,  but  used  different 

calculations  and  independent  smoothing  of  the  elastic  peaks 

to  determine  G   .   The  values  of  G    determined  by  each  of 
n  n 
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these  methods  showed  very  good  agreement,  well  within  the 

statistical  errors.   This  agreement  is  convincing  evidence 

that  the  computer  programs  used  in  the  data  reduction  did  not 

contain  errors  which  could  bias  the  conclusions  reported  in 

this  paper.   The  ratio  R  reported  in  Table  V-I  is  the  average 

R  for  those  15  runs. 

The  average  value  of  G    at  each  q  ,  using  Lomon-Feshbach 

n 
Model  #1,  is  given  in  Table  V-I I. 

2 
The  average  value  of  G    at  each  q  ,  using  Lomon-Feshbach 

n 

Model  #5,  is  given  in  Table  V-III. 

The  value  of  G    calculated  from  L-F  Model  15  and  its 

n 

S~    for  each  of  the  21  runs  is  given  in  Table  V-IV. 

GE 
n 

2 

The  average  value  of  G_,   at  each  q  ,  using  Lomon-Feshbach 

n 
Model  #15  is  given  in  Table  V-V. 

B.  ROSENBLUTH  PLOT 

2 
At  q"  -  0.2  deuteron  and  hydrogen  data  were  taken  at  60, 

7  5,  90,  and  120  degrees  in  order  that  a  Rosenbluth  Plot  could 

be  constructed.   The  theoretical  interpretation  of  the 

Rosenbluth  Plot  is  discussed  in  Section  II-G.   Table  V-VI 

gives  the  experimental  data  and  the  derived  quantities  used 

as  data  for  the  Rosenbluth  plot  given  in  Figure  V-5. 

C.  CHARGE    RADIUS    OF    THE    DEUTERON 

The  deuteron  charge  rms  radius  was  calculated  by  Method 
II  given  by  Schumacher  and  Bethe  [44] .   They  have  provided 
the  charge  radius  of  the  deuteron  from  several  deuteron  wave 
functions.   L-F  #1  and  L-F  #5  have  a  charge  radius  of  1.971  fm 
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TABLE  V-II 

2 
Average  GE   at  Each  q *  Using  Lomon-Feshbach  Model  #1 

n 


2 

q-  2 
(fm  Z) 

CE(LF#1) 

GE 
n 

S^   (Total 

\ 

0.05 

0.96881 

0.0002 

0.0015 

0.10 

0.93967 

0.0036 

0.0011 

0.20 

0.88667 

0.0028 

0.0011 

0.25 

0.86247 

0.0070 

0.0027 

0.30 

0.83961 

0.007  3 

0.0014 

0.35 

0.81796 

0.0095 

0.0019 

A  linear  fit  to  the  data  in  Table  V-II  gives: 

GE  =  (-0.0003  +  0.0012)  +  (0.0245  +  0.0059)q2,       (V-l) 
n 

while  a  linear  fit  of  G„,   weighted  by  the  statistical  S„ 

Ln  GE 

,    •  n 

only  gives: 

G^  =  (0.0001  +  0.0010)  +  (0.0231  +  0.0049)q2,        (V-2) 
E 

n 
and  the  21  points,  if  equally  weighted,  give 

G   =  (-0.0004  +  0.0016)  +  (0.0263  +  0.0072)q2        (V-3) 
E 
n 

A  plot  of  the  data  in  Table  V-II  is  given  in  Figure  V-l. 
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Figure  V-I 
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TABLE  V-III 

2 

Average  G    at  Each  q   using  Lomon-Feshbach  Model  #5 


q2 
(fm  2) 

CE(LF#5.) 

GE 

n 

S_   (Total 

GE 
n 

0.05 

0.96881 

0.0002 

0.0015 

0.10 

0.93966 

0.0036 

0.0011 

0.20 

0.88665 

0.0028 

0.0011 

0.25 

0.86244 

0.0070 

0.0027 

0.30 

0.83956 

0.0074 

0.0014 

0.35 

0.81790 

0.0096 

0.0019 

n 


A  linear  fit  to  the  data  in  Table  V-III  gives: 

G   =  (-0.0003  +  0.0012)  +  (0.0249  +  0.0059)q2        (V-4) 
En 

While  a  linear  fit  of  GL,   weighted  by  the  statistical  S„ 

n  E. 

only  gives: 

G^  =  (0.0001  +  0.0010)  +  (0.0234  +  0.0049)q2         (V-5) 
En 

and  the  21  points,  if  equally  weighted,  give 

G   =  (-0.0003  +  0.0016)  +  (0.0266  +  0.0072)q2.        (V-6) 
En 

A  plot  of  the  data  in  Table  V-III  is  given  in  Figure  V-2, 
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Figure  V-2 


106 


and  S. 


TABLE  V-IV 
(Total)  for  Each  of  the  21  Runs 


n 


(fm  2) 


0.05 


0.10 


0.20 


0.25 

0.30 

0.35 


Calculated  from  Lomon-Feshbach  Model  #15 


CE(LF#15) 


0.96899 


0.94002 


0.88731 


0.86324 

0.84049 

0.81895 


n 

0.0003 
-0.0005 

0.0079 
0.0047 
0.0011 
-0.0012 

0.0031 
0.0159 
0.0025 
0.0053 

-0.0002 
0.0023 

-0.0026 
0.0015 

0.0061 

0.0152 
0.0031 
0.0050 
0.0024 

0.0062 
0-.0107 


S_   (Total) 

X 

0.0019 
0.0023 


0.0021 
0.0017 
0.0024 
0.0021 

0.0058 
0.0051 
0.0040 
0.0024 
0.0024 
0.0025 
0.0025 
0.0037 

0.0031 

0.0028 
0.0029 
0.0034 
0.0024 

0.0026 
0.0027 


A  plot  of  the  data  in  Table  V-IV  is  given  in  Figure  V-3. 
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Figure  V-3 
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TABLE  V~V 


Average  G    at  Each  q  Using  Lomon-Feshbach  Model  #15 


2 

q 

(fm"2) 


0.05 
0.10 
0.20 
0.25 
0.30 
0.35 


CE(LF#15) 


0.96899 
0.94002 
0.88731 
0.86324 
0.84049 
0.81895 


n 

0.0000 
0.0034 
0.0021 
0.0061 
0.0062 
0.0084 


Sr       (Total) 

\ 

0.0015 
0.0010 
0.0011 
0.0031 
0.0014 
0.0019 


A  linear  fit  to  the  data  in  Table  V-V  gives: 
G   =  (-0.0000  +  0.0012)  +  (0.0200  +  0.0058)q2 


(V-7) 


While  a  linear  fit  of  G   weighted  by  the  statistical  S 

n  ( 

only  gives: 


G^  =  (0.0000  +  0.0010)  +  (0.0194  +  0.0050)q 
En 

and  the  21  points ,    if  equally  weighted,  give: 
G^  =  (-0.0003  +  0.0016)  +  (0.0231  +  0.0072)q2 


(V-8) 


(V-9) 


A  plot  of  the  data  in  Table  V-V  is  given  in  Figure  V-4, 
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of 

2 
Deuteron  Form  Factor  Squared  versus  1+2(1+7))  tan  (9/2 
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Figure  V-5 
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while  L-F  #15  has  r   =  1.965  fm.   The  experimental  ---„-  G„  (0) 

dq     n 
which  were  determined  using  L-F  #1,  5,  and  15  are  plotted 

2 

against  the  square  of  the  deuteron  charge  radius,  r„  ,  as 

bd 
given  by  Schumacher  and  Bethe,  in  Figure  V-6.   Schumacher 

notes  that  such  a  plot  yields  an  excellent  straight  line 
for  all  models  with  the  correct  binding  energy.   The  inter- 
section of  this  line  with  the  value  of ^    G„  (0)  determined 

dq     n 
by  the  very  accurate  neutron-electron  scattering  work  yields 

2 

r„  .   The  value  of  the  deuteron  charge  rms  radius  was  deter- 

Ed 
mined  from  Figure  V-6  to  be  r^   =  1.9641  +  0.0074  fm. 

Ed 
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VI.   CONCLUSIONS 


From  the  75  measurements  of  the  ratio  of  the  elastic 

electron-deuteron  cross  section  to  the  elastic  electron-proton 

2 
cross  section  at  lov;  q   the  following  conclusions  are  drawn: 

i)   the  discrepancy  in  the  neutron  charge  form  factor 
measured  in  several  earlier  experiments  [1,2,3,4]  was  re- 
solved.  The  value  of  t?   G„  (0)  was  determined  to  be  in 

dq  .    n 
agreement  with  the  thermal  neutron  results  for  all  three 

deuteron  wave  functions  (L-F  1,5,15)  which  were  used  to  analyze 

the  data. 

ii)   the  value  of  — ~   G   (0)  ranged  from  (0.0245  + 

2  dq     n 

0.0059)  fm  when  L-F  1  was  used  to  extract  G„   from  the 

hi 

2  n 

deuteron  data,  through  (0.0249  +  0.0059)  fm  when  L-F  5  was 

2 
used,  to  (0.0200  +  0.0058)  fm  when  L-F  15  was  used.   Thus 

L-F  15  produced  a  value  of  -  G   (0)  in  best  agreement  with 

dq     n  d 

the  very  accurate  thermal  neutron  results  of  — ^   Gp  (0) =0.0193+ 

2  dq     n 

0.0004  fm   and  was  determined  to  be  the  model  (of  the  three 

considered)  which  best  describes  the  deuteron.   The  method  of 
calculating  the  LF  wave  functions  for  different  models  and 
the  more  significant  differences  between  the  three  models  con- 
sidered here  is  discussed  in  Appendix  C. 

iii)   the  charge  radius  of  the  deuteron  was  measured  to 

be  r„  =  1.9641  +  0.0074  fm.   Prior  to  this  result  the  most 

Ed 
recent  measurement  [21]  was  r   =  1.95  +  0.02  fm.   The  very 

Ed 
accurate  value  of  the  reported  charge  radius  of  the  deuteron 

was  possible  because  all  data  were  taken  at  or  below 
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2  -2 

q   =0.35  fm    and  thus  it  could  all  be  used  to  determine  the 

2 

intercept  and  coefficient  of  the  q   terra  in  the  equation 

2 

G       (q2)  =  G   (0)  -  qr        r2 
E   A       E        6     rms 
n  n 

iv)  as  a  check  on  the  results  of  this  experiment  a 

2  -2 

Rosenbluth  plot  of  all  of  the  data  taken  at  q   =0.2  fm 

was  produced  in  Figure  V-5.   The  slope  of  this  line  is 

(0.0026  +  0.0012).   The  theoretical  slope  of  this  line  is 

2    2 

4rTj  G-,  (see  Section  JI-G)  .   Using  the  "scaling  law"  (Equation 

2        M    2 
11-31)  this  slope  becomes  jT?  (4J3  rp  )   where  GQ=  (G   +G   )C  . 

2  P  P    n 

At  q   =  0.2  the  theoretical  slope  is  equal  to  0.000822.   The 

experimental  slope  is  about  1.5  0_         .     from  this  value. 
1  *  Total 

However,  Table  V-VI  shoves,  as  does  Figure  V-5,  that  the 
second  120   point  has  a  significantly  different  value  of  G, 
than  do  the  others.   If  this  one  point  is  dropped,  a  linear 
fit  to  the  remaining  7  points  gives  a  value  of  the  slope  of 
0.0013  +  0.0013  in  good  agreement  with  the  theoretical  slope. 

This  shows  that  the  one  photon  exchange  is  valid  for  the 

2 

deuteron  at  this  q  . 
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APPENDIX  A 


DE  VRIES  B*  FIT 


De  Vries  et_  al_.  [10]  analysed  the  proton  cross  section 
measurements  of  BumiHer  et_  al .  [11]  which  consisted  of  58 
data  points,  the  114  points  of  Janssens  et_  al.  [10],  and  the 
6  points  of  Chen  et_  al.  [12]  as  well  as  the  71  deuteron  cross 
section  measurements  which  de  Vries  presented.   These  data 
were  used  to  obtain  the  best  fit  to  the  three  pole  Clementel- 
Villi  [40]  model  for  the  isotopic  form  factors  given  in 

Equations  A-l.   This  de  Vries  fit  is  valid  only  in  the  low 

2 

q   region. 

SE  SE 

Gg  =  0.50  {   12  +    2  2        +  (I"  S   -  S   )} 

E  1+  qZ 1+  q_ El    E2  J 

15.6  26.6 

G^  =  0.50  { K      +     (1+V-  )} 

1+  si  1 

4 

15.6       26.6 

VM 
GM=  2-353  {j^2    +  U  "VMl^  tA_1) 

4 
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The  proton  form  factors  are 

p 

GK„  =  G7„  +  gX 

M     M     M 
P 

while  the  neutron  form  factors  are 

GE  "  GE  "  GE 
n 

Gx.  =  G.,  -  G._ 

M     M     M  . 
n 

De  Vries  started  with  this  theoretical  model  given  by 
Equations  A-l  and  computed  cross  sections  which  were  then 
compared  with  the  experimental  data.   The  free  parameters 

appearing  in  the  model  were  adjusted  so  that  a  minimum  in 

2 

X   was  obtained. 

The  mass  of  the  p  meson  was  taken  as  a  free  parameter 
since  the  width  of  the  resonance  was  so  wide  (100  MeV) ,  but 
the  seven  other  free  parameters  were  reduced  to  six  by  ap- 
plying the  constraint 

■2 


(-7T  ge  )  - 


'E  /  2   ' 

dq     n  q  =0 

Data  from  the  several  experiments  mentioned  above  were 
combined  in  several  different  combinations.   A  reduced  chi- 
square  was  calculated  for  each  of  these.  'Combinations  b  and 
b'  were  recommended  by  the  authors  as  the  best  fits.   In  this 
experiment  the  values  of  the  free  parameters  given  by  the  b* 
fit  were  used  since  Buchanan  [68]  has  reanalyzed  all  known 
proton,  neutron,  and  deuteron  data  and  concluded  that 
de  Vries  b*  fit  best  describes  the  proton  form  factors.   The 
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b'  combination  includes  the  data  of  '  Janssens  et  aJL . ,  de  Vries 
et  al. ,  and  Chen  et  al. ,  and  gives  the  following  for  the 
values  of  the  parameters. 


SE  =  2.628  SM  = 


1  1 


S^   =-1.853  SM  =-3.435 

E2  M2 


V^  =  1.191  V„  =  1.064 

El  Ml 


M  2=  8.463 

r 
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APPENDIX  B 
THE  VI RIAL  COEFFICIENTS  OF  HYDROGEN  AND  DEUTERIUM 

In  this  experiment  the  target  gases  were  contained  in  a 
metal  gas  chamber  that  was  in  good  thermal  contact  with  a 
liquid  nitrogen  reservoir  which  was  at  atmospheric  pressure. 
The  chamber  was  filled  from  high  pressure  bottles  of  H?  or 
D-  through  reducing  and  regulating  valves  to  a  gauge  pressure 
of  about  150  psig.   The  density  of  the  gas  in  the  target  was 
then  computed  using  the  virial  coefficients  of  H~  and  D,-,  at 
77.35  K  (Liquid  nitrogen  temperature).   Since  the  pressure  was 
monitored  continuously  throughout  a  data  run  the  virial  ex- 
pansion was  calculated  in  terms  of  pressure  coefficients. 

Using  the  data  of  Michels  and  Goudeket  [69]  and  the  method 
given  in  the  Handbook  of  Physics  [70]  the  gas  density  was  de- 
termined to  be 

PM  f        1        >    ,      3 
RT  \+    B'P  +  C  P2 

where 

B'  =  -  0.189409  x  10~3  in2/lb 

C  =    0.1298969  x  10~7  in4/±b2  . 

3 

At  150  psig  hydrogen  had  a  density  of  0.0036723  gm/cm 

and  deuterium  0.007  3391  gm/cm  . 

During  a  run  the  pressure  was  recorded  and  was  then  used 
to  calculate  (in  FORTRAN  program  TOPRADCR)  the  gas  density 
which  in  turn  was  used  in  Equation  IV-2  to  obtain  the  number 

120 


2 
of  target  atoms  per  cm   and  in  the  bremsstrahlung  and  Landau 

radiation  correction  calculations  as  well  as  to  partially 

determine  the  cross  section  in  Equation  IV--13. 

The  American  Institute  of  Physics  Handbook,  Third  Edition, 

which  was  published  coincidently  with  the  completion  of  the 

writing  of  this  thesis,  has  virial  coefficients  different 

from  those  given  in  Reference  70.   The  data  of  Reference  70 

was  used  for  all  calculations  in  this  work.   If  the  data  of 

the  AIPH,  Third  Edition,  are  used  the  primary  results  of  this 

work  will  not  be  changed.   That  is,  the  slope  of  G   ,  the 

n 
radius  of  the  deuteron  and  the  slope  of  the  Rosenbluth  plot 

are  unchanged.   There  will,  however,  be  a  difference  in  the 

2 

intercepts  of  both  G   (0)  versus  q   and  of  the  Rosenbluth 

n 
plot.   This  change  in  virial  coefficients  will  cause  the 

intercept  of  G    to  be  raised  to  0.002  3,  which  might  in- 

n 
dicate  that  the  "contaminent"  in  the  H„  discussed  in  Section 

IV-D  is  several  tenths  of  a  percent. 
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APPENDIX  C 

LOMON-FESHBACH  DEUTERON  ELECTRIC  CHARGE  STRUCTURE  FACTORS 

The  deutcron  electric  charge  structure  factor,  Cg,  is 
defined  as 


CE  =  $[u2(r)  +  v2(r)>0  (f   )  dr 


where  u(r)  and  v(r)  are  respectively  the  deuteron  S-and  D- 
state  wave  functions  normalized  so  that  C   =  1  when  q  =  0, 
J   is  the  spherical  Bessel  function  of  order  zero,  q  is  the 
four -momentum  transferred  and  r  is  the  separation  between 
the  neutron  and  proton.   In  order  that  one  might  extract  the 
neutron  form  factor  from  data  which  give  the  shape  of  the 
deuteron  one  needs  to  have  C^,  which  means  that  one  must 
have  u(r)  and  v(r).   These  are  determined  from  the  nuclear 
model  used.   Lomon  and  Feshbach  derived  nuclear  models  from 
nucleon-nucleon  data  by  fitting  a  boundary  condition  model 
interaction  determined  largely  by  field-- theoretic  forms  [15]. 
Lomon  and  Feshbach  use  an  energy  independent  boundary  condition 
for  the  deuteron  wave  function, 

where  r   is  the  range  of  the  potential.   This  condition  must 

be  different  in  each  angular  momentum  state  [72] .   E.L.  Lomon 

[14]  has  provided  values  of  C   based  upon  three  different 

E 

nuclear  models.   These  values  of  Cp  are  listed  in  Table  C-I . 

E 
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The  three  models  which  have  been  considered  are  called  LBn#l, 
LF#5  and  LF#15  in  agreement  with  Ref.  15.   Model  L-F  #1  has 
a  deuteron  binding  energy  of  2.2242  MeV  and  4.57  percent 
D-state;  model  L-F  #5  has  a  binding  energy  of  2.2245  MeV  and 
5.202  percent  D-state;  model  L-F  #15  has  a  binding  energy  of 
2.2236  MeV  and  7.55  percent  D-state. 

TABLE  C-I 

LOMON-FESHBACH  DEUTERON  ELECTRIC 
CHARGE  STRUCTURE  FACTORS 


qV-2) 

CE(LF#1) 

CE(LF#5) 

CE(LF#15) 

0.025 

0.98414 

0.98413 

0.98423 

0.05 

0.96881 

0.96881 

0.96899 

0.10 

0.93967 

0.93966 

0.94002 

0.15 

0.91235 

0.91234 

0.91285 

0.20 

0.88667 

0.88665 

0.88731 

0.25 

0.86247 

0.86244 

0.86324 

0.30 

0.83961 

0.83956 

0.84049 

0.35 

0.81796 

0.81790 

0.81895 

0.40 

0.79741 

0.79735 

0.79851 

0.45 

0.77789 

0.77782 

0.77907 

0.50 

0.75929 

0.75921 

0.76057 

0.55 

0.74156 

0.74147 

0.74291 

0.60 

0.72462 

0.72452 

0.72604 

1.3 

0.54649 

0.54626 

0.54854 

2.0 

0.43328 

0.43293 

0.43555 

3.0 

0.32637 

0.32588 

0.32867 

4.0 

0.25438 

0.25378 

0.25657 

123 


APPENDIX  D 


THE  GAS  TARGET  SYSTEM 


A  schematic  diagram  of  the  gas  target  system  as  it  was 
developed  for  this  experiment  is  shown  in  Figure  D-l.   A 
detailed  scale  drawing  of  the  gas  target  assembly  is  given  by 
Savage  [71] .   In  this  section  both  the  best  and  the  poorest 
techniques  which  were  use'd  in  this  experiment  to  operate  the 
gas  target  system  are  discussed  and  hopefully  this  will  pre- 
vent others  from  following  some  unprofitable  techniques. 

As  the  linear  accelerator  was  being  turned  on  for  a  run 
the  gas  fill  system  was  connected  together  between  valves 
V7  and  V8  and  valves  VI,  V2  and  V7  were  opened.   (Actually 
V7  and  V8  were  always  open  during  all  the  runs  done  with  the 
final  technique.)   The  pressure  gauge  was  checked  to  see  that 
it  read  zero  gauge  pressure.   Then  valve  VI  was  closed  and 
valves  V3,  V5  and  V8  were  opened.   The  forepump  was  turned  on 
to  evacuate  the  complete  gas  fill  system  up  to  valves  V4  and 
V6  which  are  parts  of  the  gas  bottles.   This  pulled  the  system 
down  to  about  50  microns. 

After  the  accelerator  had  been  tuned  for  the  correct 
energy  a  large  (50  liter)  liquid  nitrogen  dewar  was  connected 
to  feed  the  liquid  nitrogen  reservoir.   Temperature  sensors 
were  placed  in  this  reservoir.   These  sensors  were  taped  to  a 
plastic  rod  and  were  arranged  along  this  rod  as  shown  in 
Figure  D-l.   Whenever  the  nitrogen   level  fell  below  the  top 
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Figure   D-I 
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sensor  the  Cryomiser  would  open  the  valve  in  the  compressed 
air  line  forcing  air  into  the  dewar  and,  since  the  dewar  had 
a  pressure  seal  at  the  top,  this  forced  liquid  nitrogen  into 
the  gas  target  reservoir.   Through  trial  and  error  it  was 
found  that  one  minute  of  pumping  was  satisfactory;  if  longer 
"fill"  times  were  set  on  the  Cryomiser  it  was  found  that  the 
gas  target  assembly  became  so  cold  that  it  would  harden  the 
"0"  ring  seal  and  a  serious  vacuum  leak  would  occur.   A 
shorter  fill  time  did  not  adequately  fill  the  reservoir. 
The  lower  sensor  would  sound  an  alarm  in  the  end  station  if 
the  nitrogen  level  fell  below  it,  and  this  audio  alarm  was 
reproduced  in  the  control  room. 

In  the  early  runs  the  gas  target  was  allowed  to  cool  for 
at  least  15  minutes  before  the  first  gas  was  put  in  the  target, 
In  later  runs,  the  gas  was  put  in  immediately  after  the  liquid 
nitrogen  was  pumping;  both  methods  were  satisfactory. 

After  a  vacuum  of  at  least  50  microns  was  read  on  the 
thermocouple  and  the  liquid  nitrogen  had  started  to  pump, 
valve  V2  was  closed.   A  pressure  of  50  microns  insured  that 

any  contaminants  in  the  gas  fill  system  would  add  less  than 

5 
one  contaminant  molecule  in  10   target  molecules.   Both  regu- 
lators were  turned  to  the  closed  position  and  valves  V3  and 
V5  would  also  be  closed.   Valve  V4  was  then  opened.   The 
pressure  gauge  marked  HP,  for  high  pressure,  read  between 
2000  to  200  psig  depending  of  course  on  the  amount  of  gas 
left  in  the  bottle.   The  LP,  low  pressure,  gauge  should  then 
read  zero.   If  it  did  valve  V3  was  opened.   Since  the  gas 
target  had  1  mil  stainless  steel  windows  the  regulator  was 
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moved  slowly  and  steadily  in  order  that  these  windows  would 
not  be  unnecessarily  strained  as  the  pressure  increased  from 
vacuum  to  about  150  psig.   This  pressure  was  read  on  both  the 
LP  gauge  and  on  the  much  more  accurate,  larger  gauge  in 
Figure  III-2  simply  shown  with  the  word  pressure  on  its  face. 
The  regulator  maintained  the  system  very  near  the  set  pressure 
throughout  all  runs  although  there  usually  was  some  pressure 
drop  for  a  short  time  while  the  gas  and  nitrogen  systems 
reached  equilibrium.   A  remote  control  TV  camera  was  panned 
from  the  gas  target  to  the  large  pressure  gauge  after  the 
operator  had  insured  that  the  beam  was  properly  aligned 
through  the  gas  target.   The  camera  then  generally  remained 
on  the  large  pressure  gauge  throughout  the  run.   Since  the  gas 
pressure  would  drop  whenever  a  nitrogen  fill  was  progressing, 
the  operator  recorded  the  pressure  of  the  gas  at  the  start 
of  the  fill  (i.e.,  when  the  audio  warning  went  on  in  the 
control  room)  and  he  would  also  record  the  pressure  at  the 
end  of  the  1  minute  fill  (when  the  audio  stopped) . 

After  all  the  hydrogen  data  were  taken  valve  V4  was 
closed  and  valves  V2  and  V5  were  opened.   Valve  V2  was  opened 
very  slowly  in  order  that  the  1  mil  stainless  window  would  not 
be  unloaded  too  quickly  as  it  experienced  the  decrease  from 
150  psig.  to  vacuum.   Valve  V5  was  opened -at  this  time  to 
ensure  that  no  hydrogen  leaked  through  it  during  the  hydrogen 
fill.   Again  a  vacuum  of  at  least  50  microns  was  pulled  on  the 
entire  system  up  to  valves  V4  and  V6  (the  regulators  did  not 
provide  a  vacuum  barrier) .   After  the  system  was  down  to  at 
least  50  microns,  valves  V2,  V3  and  V5  were  closed.   The 
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regulator  on  the  deuterium  was  checked  to  ensure  that  it 
was  closed,  thus  there  was  no  step-like  rise  in  pressure  on 
the  gas  target  windows  when  valve  V6  was  opened.   Again  if 
the  deuterium  LP  gauge  read  zero,  valve  V5  was  opened  and  the 
deuterium  regulator  was  turned  slowly  to  reduce  the  strain  on 
the  1  mil  steel  window.   After  reading  about  150  psig.  on 
the  large  gauge  the  system  had  some  time  to  come  to  equili- 
brium as  the  operator  re-tuned  the  accelerator.   After  checking 
the  position  and  appearance  of  the  beam  on  the  target  the  re- 
mote TV  camera  was  again  put  on  the  large  pressure  gauge  and 
the  gas  pressure  was  again  recorded  at  the  start  and  finish 
of  each  nitrogen  fill  as  the  deuterium  data  were  taken. 

After  taking  the  deuterium  data  valve  V6  was  closed  and 
valves  V2  and  V3  were  opened  and  the  system  was  again  pulled 
to  at  least  a  50  micron  vacuum.   Then  the  accelerator  beam 
was  sent  through  the  empty  target  in  order  to  obtain  empty 
target  background  at  several  spectrometer  settings  over  the 
range  of  data  taken  on  that  run. 

After  the  accelerator  was  secured  the  forepump  and  the 
cryomiser  were  also  turned  off  and  valve  VI  was  opened.   This 
let  the  gas  system  down  to  air  and  condensed  some  gases,  e.g. 
water  vapor  in  the  gas  target.   The  gas  target  was  not  left 
at  vacuum  because  it  was  believed  that  the  1  mil  windows 
would  fail  if  the  target  chamber  should  have  had  a  vacuum 
failure.   Very  early  in  these  experiments  valve  V8  was  then 
closed  in  order  that  excess  water  would  not  condense  in  the 
gas  target,  but  as  the  liquid  nitrogen  evaporated  during  the 
night  the  1  mil  windows  were  blown  out  -  presumably   due  to 
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high  vapor  pressure  within  the  gas  target.   The  later 
technique  was  to  disassemble  the  fill  line  at  valve  V8 .   This 
protected  the  target  from  a  mishap  as  described  above  should 
someone  have  accidently  closed  V8  or  even  VI. 

Finally  the  compressed  air  line  was  separated  and  was 
vented  directly  into  the  liquid  nitrogen  reservoir  directly 
above  the  gas  target.   This  accelerated  the  rate  of  evapora- 
tion of  the  liquid  nitrogen  and  generally  would  have  all  the 
nitrogen,  and  more  importantly,  all  the  water  which  condensed 
during  this  evaporation  removed  from  the  reservoir  within 
12  hours.   It  was  found  that  if  the  experimenter  starts  to 
fill  the  reservoir  with  liquid  nitrogen  for  another  run  before 
the  condensed  water  from  the  last  run  is  completely  evaporated, 
then  he  will  most  probably  cause  the  indium  seal  between  the 
gas  target  and  the  reservoir  to  fail  severely  -  in  fact  in 
two  experiments  the  target  chamber  went  to  air  and  only  the 
quick  closing  automatic  valve  at  the  collimator  saved  the 
electron  gun. 

It  should  be  mentioned  that  the  gas  target  with  valve  V8 
closed  will  not  leak  any  detectable  amount  of  gas  during  at 
least  a  10  hour  period.   In  the  early  part  of  this  work  valve 
V8  was  closed  after  the  initial  filling  of  the  target  with 
gas,  but  the  pressure  at  the  end  of  the  run  generally  did 
not  reproduce  the  initial  pressure  -  it  was  generally  higher. 
These  pressure  variations  were  of  the  order  of  one  percent  - 
too  much  of  a  variation  to  permit  the  accuracy  of  measurement 
required  for  these  experiments.   Thus  the  open  target  method 
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described  above  was  used  with  continuous  monitoring  of  the 
pressure  throughout  the  run,  i.e.  by  recording  the  pressure 
at  the  start  and  finish  of  each  nitrogen  fill. 


130 


APPENDIX  E 
COMPUTER  PROGRAMS  USED  IN  DATA  REDUCTION 

The  following  FORTRAN  programs  were  used  to  calculate 
quantities  explained  elsewhere  in  this  work. 

TOPRADCR  calculated  the  Schwinger,  bremsstrahlung  and 
ionization  (Landau)  radiation  corrections  as  well  as  the 
density  of  the  deuterium  and  hydrogen. 

TOP#FORM  calculated  the  experimental  cross  section  for 
each  channel.   The  sample  listed  is  the  program  for  channel 
number  5  (i.e.  r\   =5).   It  was  also  used  to  obtain  the 
smoothed  data  which  was  plotted  to  define  the  elastic  peak. 

TOPROTON  calculated  the  proton's  Mott  and  Rosenbluth 
cross  sections  of  the  proton,  the  magnetic  correction  term 
and  the  electric  and  magnetic  form  factors  of  the  proton 
based  on  de  Vries'  b1  fit. 

TOPDEUTR  calculated  the  Mott  cross  section  of  the  deuteron 
and  the  magnetic  correction  term  which  includes  within  it  the 
relativistic  correction  to  the  wave  functions. 
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